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NONANALYTIC FUNCTION GENERATION ROUTINES 
FOR 16-BIT MICROPROCESSORS 
by James F. Soeder and Maryrita Shaufl 
Lewis Research Center 


SUMMARY 

This report describes various interpolation techniques for three types 
of nonanalytic functions: univariate, bivariate, and map. These interpola- 

tion techniques are then implemented in scaled-fraction arithmetic on a 
representative 16-bit microprocessor. This work was done on an Intel 8086; 
however, the programs can be modified for use with any 16— bit microprocessor. 
A Fortran program is described that facilitates the scaling, documentation, 
and organization of data for use by these routines. Listings of all these 
programs are included in an appendix to the report. 


INTRODUCTION 

As microprocessors become more sophisticated, they will be used in in- 
creasingly complex applications. Specifically, advanced third— generation 
microprocessors such as the Zilog 8000, the Motorola 68000, and the Intel 
8086 look more like minicomputers than like the programmable digital con- 
trollers that were characteristic of their predecessors. Therefore these 
third— generation microprocessors will be called on to control and simulate 
systems that require increasingly complex calculations in less and less time. 
One common type of calculation that is required in the simulation and control 
of a system such as a gas turbine engine is high-speed nonanalytic function 
generation. A nonanalytic function is taken here to mean any function of one 
or two variables that can be described by a table of values. 

This function generation is accomplished by a program that can inter- 
Pof^te univariate, bivariate, and map functions rapidly enough to provide 
stable operation of the system. In many cases the generation of nonanalytic 
functions must be done as fast as possible. The current third— generation 
microprocessors can do this most economically by using assembly language and 
scaled-fraction arithmetic. Although the time required to write and debug 
these assembly language programs is somewhat greater than that required for 
high-level language programs, the resulting programs will be faster. In 
addition, if their calling sequences are general enough, they can be used 
in other applications, even those including high-level languages. The draw- 
back of using scaled-fraction arithmetic can be overcome very easily by using 
a Fortran program to scale data. The program can accept the function data in 
engineering unit form, scale the data, and output them to a file that can be 
directly assembled by an assembler. By processing data in this automatic 
manner manual transcription and scaling errors can be virtually eliminated. 
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This report describes three general function types: univariate, bivari- 

ate, and map. It then characterizes the interpolation techniques for each 
type of function. These interpolation techniques are similar to the techni- 
ques discussed in reference 1. Next, assembly language routines for function 
generation, which were written for the Intel 8086, are examined. Although the 
routines are specific to this 16-bit microprocessor, it is merely used as a 
demonstration vehicle. The ideas and concepts can be applied to any micro- 
processor with a 16-bit architecture and fixed-point-hardware, multiply-and- 
divide capability. Finally, the Fortran program that takes data in engineer- 
ing units and converts them to scaled-fraction form is described. Listings 
of all these programs are included in the appendix. 


DESCRIPTION OF NONANALYTIC FUNCTION TYPES 

There are three basic types of functions that are considered in advanced 
microprocessor applications: Univariate, bivariate, and map. Descriptions of 

each function and their various interpolation techniques follow. 


Univariate Function 


A univariate function is the simplest type of nonlinear relationship. 
Figure 1 shows an example of a univariate function. It merely consists of a 
set of output values y, corresponding to a set of input values x. The com- 
putation of an intermediate y value corresponding to a particular x is done 
simply by using the linear interpolation equation 




( y H ' y L> + y L 


(D 


where the definitions of the x T s and y’s are shown in figure 1. This 
simple interpolation equation can be implemented in one of two ways. First, 
one can store a set of scaled x and y points and compute equation (1) 
directly. Second, one can rewrite equation (1) in the following form: 


y v = m (xy) + b 


( 2 ) 


where 



and 


b = y intercept 

In other words, for each function segment that is defined by two breakpoints, 
a linear equation can be written. The slope of the segment is computed in a 
straightforward manner, and the y intercept is the point where the extension 
of the segment under consideration would intersect the y axis. Once the m's 
and b r s are computed off line, they are stored along with the x f s and then 
used with equation (2) to generate the function. For simplicity of reference 
the first method is called FUN1 and the second NEFG. 
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FUN1 has the advantages of being straightforward, being easy to scale, 
and using a smaller amount of memory storage than NEFG. NEFG has the advan- 
tage of being faster because the divide operation necessary to compute the 
slope has been eliminated. However, NEFG does require considerably more 
storage per function and more off-line computation to determine the correct 
m's and b's. The tedium of this computation can be eliminated by using 
the Fortran program described later, in the section FORTRAN DATA PROCESSING. 


Bivariate Function 


An example of a specialized bivariate function is presented in figure 2. 
In this nonlinear function the output value y is dependent on two inputs, 
x and z. The curves, each of which has a particular value of z, have the 
same number of breakpoints . However, breakpoints occur at the same values of 
x. This simplifies the interpolation of this bivariate function. Therefore 
one only needs to perform a FUN1 type of univariate function routine on each 
of the z curves in order to find new y values for the input x and then 
to interpolate the two y values by using the z input. This procedure can 
be summarized by the following equations and reference to figure 3: 


yv ~ f( x v ,z v) 

Y L = (xg - x") (y B " y A) + FA 
/x v - x L \ 

y n " [iq - x L J (y D - yc) + yc 

_ ( z v - z l\ , 

Yv = Vn - z lJ (y H ' Yl) + y L 


(3) 

(4) 

(5) 


where all the x's, y's, and z’s are defined in figure 3. These equations 
can be manipulated in a straightforward manner to produce the output of the 
specialized bivariate function. Implementation of this method is called 
FUN2. A method could be devised that uses an NEFG interpolation technique 
to generate y L and y H . However, the scaling and storage requirements 
become quite cumbersome. This negates the time advantage gained by elimin- 
ating the divide operation required in equations (3) and (4). 


Map Function 

An example of a map function is shown in figure 4. These functions are 
similar to the bivariate ones because they have the same number of points for 
every z curve. However, in map functions each curve can have a unique set 
of x. and y breakpoints. Because of this, one cannot use the simple inter- 
polation technique that worked for the bivariate function; that is, since the 
x-coordinate did not change from curve to curve, one only had to compute two 
new y values and interpolate on z. However, for a map function, one must 
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compute the position that a curve would have if it had the desired value of z 
This is done by implementing the following series of five equations with ref- 
erence to figure 5: 


y v = f ( x v» z v) 

X F = X C + ( zh - (x B " X C> (6) 

X C = X D + ( Z V R - z£) (x E " X D> (7) 

yp = y c + (-^4^) (y B - y c> (8) 

y G = y D + (zh ~ z ' l ) (y E “ V (9) 

y V = y F + (xg - 4) (Yg " Yp) (10) 

where all the x's, y's, and z's are as defined in figure 5. These equa- 
tions can then be implemented directly. This interpolation technique is 
referred to as FUN3. 


DESCRIPTION OF ROUTINE IMPLEMENTATION 

The assembly language implementation of the four interpolation tech- 
niques is done in a straightforward manner. These routines can be called 
either from an assembly language program or from a high-level language 
program written in Intel f s PL/M-86 (refs. 2 and 3). In the former case, 
parameters are passed back and forth by using registers. In the latter, 
parameters are passed to the routine on the stack and returned to the call- 
ing program by the AX register (the accumulator). These routines were 
written for the Intel 8086 implementations, which use up to 64K of memory. 
Some changes would be necessary if the programs are to be used in the full 
1-megabyte memory environment. Complete listings of the programs are given 
in the appendix and individual details are covered below. 


Univariate Functions (FUN1) 

The FUN1 routine is the implementation of the first method of univariate 
function lookup. Several features of the program are worth mentioning. The 
first feature is the format used to store the data, an example of which is 
shown in figure 6. By having minus and plus full scale at each end of the 
data there is no need to selectively limit the x input. However, a maximum 
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check is made to make sure the input x value is not equal to 32767 and, 
if it is, to decrease the value by 1 bit. This prevents the routine from 
getting lost looking for a value that is greater than or equal to 32767. 

(See FUN1 label point in program listing.) 

The second feature is that the routine has three entry points: PFUN1, 
FUN1, and FUN1A. PFUN1 is used to pass parameters by using the stack. This 
is the method that would be used from a call by a PL/M-86 program. FUN1 is 
used when parameters are passed through the registers, as would be done in 
an assembly language call. FTJN1A should be used as an entry point to look 
up another function that has the same x breakpoints but different y break- 
points (as shown by the example format of fig. 7), and the x input values 
for both the functions are the same. With the FUN1A call the next function 
can be looked up by using the index values and slope factors that have already 
been passed to, or computed by, the program. This is useful because, once one 
computes the interpolation factor of equation (1) , much of the computation is 
done. 


The third feature is that this routine passes a new table pointer to the 
ca lli n 8 program along with the interpolated value. By passing a new table 
pointer to the calling program, on the next call of this particular function 
the FUN1 program will not have to search from the beginning of the table for 
the upper and lower bound of the interpolation. This makes the overall inter- 
polation time less since the present interpolation interval will probably be 
within the same interval or one interval away from the last interpolation 
interval. The passing of a table's x index value is done indirectly from 
within the routine. This is the violation of the PL/M-86 compiler's block 
structure rules governing the change of formal calling parameters by an 
external procedure (ref. 3). However, if one is aware that this updating is 
done, it should pose no particular problem, and a good deal of computa- 
tion time will be saved. 

Finally, shifts are made in the program at strategic places to assure 
that when dealing with fixed-point numbers no overflows will occur and produce 
erroneous results. This can happen in two different cases. First, if the 
most significant half of the dividend is greater than, or equal to, half the 
divisor, an overflow will occur. Second, if two adjacent x breakpoints are 
present that have an absolute distance greater than 32767 (i.e., -18000 and 
17000), a subtract overflow will occur. The first overflow problem can be 
eliminated by right— arithmetic shifting the dividend an appropriate number 
of places before the divide. The second can be eliminated by noting that the 
interpolation factor must always be positive. Therefore, if the interpola- 
tion factor is computed by using an unsigned divide, subtract overflows make 
no difference. 


Univariate Functions (NEFG) 

The other univariate function routine, NEFG, has many of the same fea- 
tures as the FUN1 routine. It has a calling sequence for a PL/M-86 transfer, 
PNEFG , and one for an assembly language transfer, NEFG. In addition, it has 
the capability to pass back the table index value that points to the current 
interpolation interval in much the same manner as FUN1. Like the index in 
FUN1, this provides for speedup the next time the function is called. 
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The data storage method in this routine, however, is quite different 
from that in FUN1. Figure 8 presents an example of the KEFG data storage 
format. The x and y values are single precision. However, the b 
values are double precision and therefore take up two word locations. The 
reason for the double-precision storage can be seen from the scaling that 
must be done: 

y (!) ' m (t) x (!) + b (!) <n 

where (A/C) , (D/E), and (F/G) are the scale factors of the respective vari- 
ables. Reflecting on this relationship, since 



( 12 ) 

(13) 

(14) 

(15) 


must be true. However, it cannot always be guaranteed that conditions (14) 
and (15) can be met. Therefore the scale factors must be adjusted by a 
factor k: 

-32768 i m ^ A <. 32767 
-32768 < b (A) I < 32767 


The k factor can be made anything, but making it a power of 2 allows simple 
shifting to implement it. The exact calculation sequence for computing the 
k’s for each function is described in the section Scaling Section; but if 
the sequence yields a k other than 1, the final answer must be left-shifted 
that many times outside the routine. This will not yield an overflow since 
by definition the final answer y must lie between -32768 and 32767. There- 
fore it is necessary to make the b values double precision because sometimes 
the number of shifts becomes as large as five or six. Making this many shifts 
in single precision could cause an intolerable reduction in accuracy. In an 
assembly language call the required double-word left shifts can be implemented 
directly in the calling program. For a PL/M-86 call subroutine SHFT is pro- 
vided to do the shifts the appropriate number of times. A listing of this 
routine is given in the appendix. Finally, note that there is no economy 
associated with computing two functions in the piggyback style of FUN1 and 
FUN1A. This is because of the complexity involved in computing the new m 
and b entry points. 
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Bivariate Function (FUN2) 

The FUN2 routine is a straightforward implementation of equations (3) 
to (5). This routine, like FUN1, has three entry points - FUN2, PFUN2, and 
FUN2A. The PFUN2 call uses the stack to pass all calling parameters. How- 
ever, the FUN2 call uses the registers and the stack to pass parameters and 
thereby reduce calling overhead. Since the routine uses two FUN1 lookups, 
it implements logical shifts and unsigned divides in the same critical places 
as FUN1 to avoid overflow problems. Figure 9 presents an example of the FUN 2 
data storage format. In this format each row of y's corresponds to the 
respective x s at a particular z. The first and last rows of y data are 
the same as the row that follows or precedes them, respectively. This allows 
implicit limiting of the z input in the same manner as the x input. The 
routine also employs both an x and z index value that is passed to the 
calling program. These index values allow the routine to be close to the area 
where the function interpolation will take place on the next call. Finally, a 
FUN 2 A call sequence is provided, which allows the computation of a new y if 
the input x and z values, and hence the interpolation factors, are the 
same as those for the previous curve. 


Map Function (FUN3) 

The FUN3 routine, like FUN2, is a straightforward implementation of equa- 
tions (6) to (10). This routine, also like FUN2 , has two entry points, PFUN3 
and FUN3, for various parameter passing options. However, it passes back only 
the interpolated output and the z index pointer. Since the x index pointer 
could be different depending on the two z curve values used in the interpola- 
tion, it was decided that passing two pointers back to the calling program 
would be too cumbersome for any advantage it might provide. 

An example of the data format for FUN3 is shown in figure 10. The data 
in this figure describe a map that has four z curves with three (x,y) break- 
points for each curve. Note that in this data storage format the curves are 
limited in the x direction by putting in the x value maximum and minimum 
limits (32767 and -32768). However, for the z values the limits must be 
imposed outside the routine. Otherwise, when the z search commences, the 
z pointer would be lost if the z value were above or below the maximum or 
minimum z boundary point. This technique is used to contain the pointer 
since replicating the high and low curves as is done in FUN 2 could result in 
unpredictable overflows. Finally, since map functions tend to be complex and 
to handle, most of the time one would want to start computing from 
the beginning. Therefore no provision has been made for an entry point that 
uses interpolation factors already computed (i.e., FUN1A or FUN2A). 


FORTRAN DATA PROCESSING 

The use of Fortran programs to process the data needed for the function 
generation routines makes that job much easier and quicker. Three programs 
are used to limit, scale, and format the data. They are a main (or calling) 
program, a scaling routine, and an output routine. Since the scaling and 
output routines are general routines, any user-written main program contain- 
ing data in engineering units can call them. The data for one curve are 
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passed to the scaling routine for each call. In the scaling routine the data 
and user options (other calling arguments) are analyzed, a :d the data are 
scaled and then output in both self-documenting tabular and Intel 8086 assem- 
blable form. Even though the second output format was written for the Intel 
8086, by changing a few format statements the user can easily adapt this output 
to any microprocessor, thus enhancing the flexibility of the routine. 

The user can choose several options. Curves can have "flats" added (where 
x values of -32768 and/or 32767 are added and y values from the low or high 
end, respectively, are duplicated), or they can be extrapolated on either or 
both ends. The type of function routine (NEFG, FUN1, FUN1A, FUN2, or FUN3) 
that a curve represents can be specified. Also, the relocatable or absolute 
format for assembly code generation can now be chosen. 


Main Program Data Format 

The main program, as mentioned, contains the curve data and calls to the 
scaling routine INTSCL. The data are arranged in separate arrays by curve and 
contain seven pieces of information (fig. 11) • These are (in order) x border 
specification (XB) , x breakpoint specification (XBRKP) , the number of x f s 
per z curve, the number of z curves, x values, z values (if any), and y 
values. The x border specification indicates whether flats are to be added 
or the curve extrapolated : 0.0 for flats on both ends of the curve and no 

extrapolation, 1.0 for a flat on the low end and extrapolation of the high end, 
2.0 for a flat on the high end and extrapolation of the low end, and 3.0 for 
no flats but extrapolation of both ends. The x breakpoint specification can 
have two values: 0.0 if the x breakpoint is the same for every z curve 

(i.e., NEFG, FUN1, and FUN2 functions) or 1.0 if the breakpoints are different 
(FUN3 function). The number of x’s per z curve is specified by a real 
number, as is the number of z curves. In the latter parameter provision is 
made for either 0.0 or 1.0 to represent one z curve. 

For NEFG and FUN1 curves the x, y, and z values are in a condensed 

format. The x values are listed and then the corresponding y values. No 

z values are required for the univariate function. A specific example of the 
data input format of figure 11 for a univariate curve with three points and a 
flat on each end is shown in figure 12. The formats for functions FUN2 and 
FUN3 are a little more complex. For FUN2, x values are as before listed 
first. Next the z values are listed and then the y ? s. However, the first 

row of y's corresponds to the first y value for each z curve. For FUN3 

this same pattern is followed for z and y. The x T s in this case are 
listed in the same manner as the y f s. Figures 13 and 14 represent FUN2 and 
FUN3 data input formats for three z values and four x values. In these 
figures the first subscripts correspond to the x values and the second to 
the z values . 


Calling Sequence 

Referring to the Fortran listing, data array INFO is the first parameter 
in the call to the scaling program. The other parameters, listed in order by 
their variable names, in the INTSCL subroutine are TITLE, XNUM, XDNOM, YNUM, 
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YDNOM, ZNUM, ZDNOM, FUNC, FRMT, and ORG. TITLE, the eight-character title 
of the curve, is formulated in the main program. XNUM and XDNOM, the numera- 
tor and denominator of the scale factor for the x values, are specified as 
if machine units were being converted to engineering units. That is, XNUM = 
150.0 and XDNOM = 32000.0 will give the correct scale factor (213.333) for the 
conversion of -153.600 engineering units into -32768 machine units. YNUM, 

YDNOM and ANUM, ZDNOM, the y and z scaling factors, respectively, are sim- 
ilarly specified. FUNC indicates which function routine the data array belongs 
to: 0 for NEFG, 1 for FUN1, 11 for FUN1A, 2 for FUN2 , and 3 for FUN3 . FRMT 
indicates the format of the assemblable output dataset. A zero represents a 
relocatable format and a 1, an absolute format. If 1 is specified, a value 
for the parameter ORG must also be specified, where ORG is the hexadecimal 
starting address of the data. 


Scaling Subroutine INTSCL 

When control is passed from the main program to the subroutine INTSCL, 
the actual scaling computation begins. Here the options specified by the 
user in the calling sequence are analyzed, data are scaled, and scaled data 
are output in two forms. The routine can be examined by looking at the parts 
of the program that perform the three main functions (analysis, scaling, and 
output) as three discrete sections. A functional flowchart of this program 
is shown in figure 15. 


Analysis of Calling Parameters 

The first of these program sections, the analysis section, begins by 
forming the data name arrays used on output. It then breaks the input data 
array into its several components and rearranges the x and y data. This 
rearrangement is done so that the type of curve being processed does not 
affect the routine’s treatment of the data. Finally, flats are added and 
extrapolations performed as necessary. (For FUN2 only, z automatically has 
flats added to limit the z curve values.) 


Scaling Section 

In the scaling section, scale factors are computed and z, x, and y 
values scaled. The unsealed values are saved for later use in calculation 
and output, and all scaled values are checked against the minimum and maximum 
values of —32768 and 32767. No attempt is made to adjust the scale factors, 
but out-of-range values are set to either -32768 or 32767 and the appropriate 
warning message is written to notify the user. 

The slope m and the y-intercept b values for NEFG are also calculated 
and scaled in this section. The maximum and minimum scaled values for both m 
and b are tested against the minimum and maximum of -32768 and 32767, but 
here an attempt is made to adjust the calculated scale factors by shifting and 
retesting. Thus the shift factor k mentioned in conjunction with the NEFG 
routine itself is produced. If right-shifting the original values by eight 
(dividing by 256) still does not bring the values into range, the attempt is 
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considered complete and a warning message is written to the user. Next, b 
values are separated into their most significant bits (MSB N and least signif- 
icant bits (LSB) by taking their integer and fractional portions and scaling 
them separately. The integer part multiplied by the b scale factor is con- 
sidered the most significant portion, and the fractional part multiplied by 
65536 is considered the least significant portion. 


Output Section 

The output portion of the program writes the data out in two different 
formats. The first, self -documenting, format is a series of tables that 
include the curve type and name; any warning messages; scaling factors (both 
computed values and components); and x, y, m, and b values in both engineer- 
ing units and machine units. Examples of FUN1 and NEFG tables are shown in 
figure 16. In figure 17 only the first page of the FUN 2 and FUN3 tables is 
given. These require an additional page for each z value. These tables 
give at a glance all the information necessary to characterize a curve and to 
provide documentation for permanent storage. 

The second format, shown in figure 18, is the scaled-data output in either 
relocatable or absolute format. (The absolute format differs from the relocat- 
able format only in that an ORG statement with the specified hexadecimal address 
is the first line output for a curve.) The data name - consisting of an x, y, 
z, m, or b followed by the first six characters of the curve name (five 
for FUN1A curves) and a blank - is written on the first line, along with a "DW M 
and corresponding data. The next line written has a n DW u followed by more data, 
a format that is repeated until the data for that array are exhausted. The data 
output for each of the curves is written to the dataset in the particular 
formats discussed earlier. 


Subroutine WRTDTA 

Subroutine WRTDTA is called to write all the x and y arrays for the 
Intel 8086 assemblable format. Its calling sequence has five parameters: 
FUNARY, XPTS , DNAME, FNTN, and ZPTS. FUNARY is the array of scaled data to 
be written, XPTS the number of x points in the curve, and DNAME the data 
name associated with the array. FNTN is the parameter that differentiates 
function types (0 indicates NEFG, FUN1, or FUN1A; 1 indicates FUN2 or FUN3) 
and is used primarily as a check for the number of z points, specified by 
ZPTS. 


Special-Purpose Routine 

One other routine is called from both the main program and INTSCL. This 
routine, F4MVC, is a character manipulation routine resident on the IBM 360. 

Its calling sequence parameters are SSTRNG, LI, DSTRNG, L2, and NBYTES. SSTRNG 
is the source string, LI the index of the first byte to be copied, DSTRNG the 
destination string, L2 the first byte to be replaced, and NBYTE the number of 
bytes to be copied. This routine simply moves characters, byte by byte, from 
one string to another and can therefore be replaced by a user-written Fortran 
program. 
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CONCLUDING REMARKS 

This report has described several common nonanalytic function types and 
interpolation techniques that can be used on them. In addition, a fixed-point 
arithmetic, assembly language implementation of the routines was demonstrated 
by using the Intel 8086 microprocessor. The fixed-point arithmetic was used 
for the speed and hardware minimization it currently provides. Furthermore 
the routines were written in such a way as to provide high-speed lookup while 
freeing the user from being concerned by scaling overflows. Although these 
routines were written for the Intel 8086, the overflow and scaling techniques 
are general enough concepts to be used with any 16-bit microprocessor . Fin- 
ally, a Fortran program was discussed that can be used in conjunction with 
these assembly language routines. This program can take data that are cur- 
rently in a simulation deck or any engineering unit form and scale them into 
the proper integer values that can be used by the microprocessor. The program 
then scales all the necessary function data and outputs them to a dataset that 
can be processed by the assembler. In addition, the program outputs listings 
to document all the scaled data. Therefore by using these routines one can 
realize error-free translation of function curves from engineering unit data 
to scaled-fraction data that can be used by the microprocessor for real-time 
control or simulation. 
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APPENDIX - NONANALYTIC FUNCTION GENERATION ROUTINES 
AND FORTRAN SCALING PROGRAM 



Kcs-si macro assfhps fr funi 
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3 SEPT a 


ISIS-II MCS-S6 MACRO ASSEMBLER V2.1 ASSEMBLY OF MODULE FUN1 
OBJECT MODULE PLACED IN : F ! : FFUN1 . OBJ 

ASSEMBLER INVOKED BY: ASM86 :F1:PF0M1.SRC XREF DATE (3 SEPT SO) 


IOC 


LINE SOURCE 

I ; 

•*> a 

3 ! SUBROUTINE FUN) AND FUN! A FOR UNIVARIATE 

4 ; FUNCTIONS LOOK-IJP OF NON-EQUAL INCREMENT FUNCTIONS 

5 ; 

6 ! CALL SEQUENCES FUN! 

7 ; 

8 ! SI REGISTER CONTAINS THE NUMBER OF 

9 ; POINTS IN THE TABLE 

10 5 BX REGISTER CONTAINS THE X TABLE BASE ADDRESS 

II ; DI REGISTER CONTAINS THE X TABLE INDEX VALUE ADDRESS 

12 ) AX REGISTER CONTAINS THE X VALUE 

13 ; CALL FUN1 

14 ; 

15 ; 

16 ; 

17 ? FOR FUN i A CALL DO NOT DESTROY f!X,BP,SJ,DI REGISTERS 

IS ; 

19 ; THE Y VALLE Cf THE CURVE IS RETURNED IN THE AX REGISTER 

20 ; 

21 ; IHHiHHHHH BEWARE IHHHHHIlHHtmHHHHHIHI 

t.2. * 

23 ; THE X TABLE INDEX LOCATION IS AUTOMATICALLY UPDATED BY 

24 I BY THE ROUTINE AND THEREFORE DOESN'T FOLLOW THE PL/M-86 

25 ) BLOCK STRUCTURE 

26 5 

27 : 4444444*444»444444444444444444«44iHfH444m449**4444«4*94*4 4 

23 ; 

2-9 ; CALLNG SEQUENCE FOR PL/M-86 CALL 

30 5 RSLT = F'FUNl ( X » XBSi X INDX . NOPTS ) 

31 ; 

32 ; X= VALUE OF FUNCTION 

33 ! XBS= X TABLE BASE ADDRESS 

34 ? XINDX = ADDRESS OF X TABLE INDEX VALUE 

35 ! NOPTS = NUMBER OF POINTS IN TABLE 

36 ; 

37 ; 

38 ; 

39 ; 44444444444444444444444444444449444444444444444444444 

40 ? 

41 NAME FUN1 

42 CGROUP GROUP CODE 

43 ASSUME CS: CGROUP 

44 PUBLIC FUN1,FUN1A,PFUN1 

45 ! 

46 ; 

47 ! CODE SEGMENT FOR ROUTINE 

43 ; 

49 CODE SEGMENT PUBLIC 'CODE' 

50 : PROGRAM ENTRY IF ENTERED BY PL/M-86 CALL ROUTINE 


PAGE 1 
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MCS-86 MACRO ASSEMBLER FUN1 3 SEPT 80 


LOC OBJ 

LINE S 

OURCE 



0000 5A 

51 ; 

52 PFUM15P0P 

bx ; 

TEMP SAVE OF PC 

0001 5E 

53 

POP 

SI ; 

NUMBER OF POINTS IN TABLE 

0002 5F 

54 

POP 

DI ; 

ADDRESS OF TABLE INDEX 

0003 5B 

55 

POP 

BX : 

X TABLE BASE LOCATION 

0004 58 

56 

POP 

AX : 

X TABLE VALUE TO BE LOOK-UP 

0005 52 

57 

58 

PUGH 

DX ; 

RESTORE PC TO THE STACK 


59 

60 
61 
62 

63 F 

PROGRAM ENTRY IF 
LANGUAGE ROUTINE 

CALLED BY ASSEMBLY 

0006 8B15 

UNI! MOO 

DX » 101 ] 

: MOVE ACTUAL INDEX VALUE TO DX 

000:3 87B7 

64 

65 ; 

66 i 

XCHG 

DX»DI 

; EXCHANGE TO ENABLE INDEXING 


BEGIN COMPARISON OF TABLE VALUES 

OOOA 3DFF7F 

67 

CMP 

AX, 32767 

; CHECK IF X VALLE EQUAL TO PLUS FULL SCALE 

GOOD 7501 

68 

JNZ 

CFUN 

! JUMP TO ROUTINE TO START CHECKING 

OOOF 48 

69 

DEC 

AX 

; DECREMENT X VALUE-MAKE IT =32766 

0010 3B01 

70 CFLW: CMP 

fiX > CPX+DI 1 ; COMPARE VALUE IN AX WITH TABLE LOCATION 


71 



! POINTED TO BY BX REGISTER 

0012 7D09 

72 

JOE 

ILOP 

; JUMP TO INTERPOLATION ROUTINE 


73 1 

! DECREMENT LOOP 


0014 83EF02 

74 CHOP: SUB 

DI,2 

; SUBTRACT 2 FROM TABLE POINTER 

0017 3B01 

75 

CMP 

AX,CBX+DI] ! COMPARE AX REGISTER VALLE TO TABLE LOCATION 

001? 7B0C 

76 

JGE 

INTA 

; JUMP TO INTERPOLATION ROUTINE 

001 B EBF7 

77 

JMP 

CHOP 

; JUMP TO DECREMENT LOOP 


78 ! 

; INCREMENT LOOP 


00 IB 83C702 

7? 

ILOP: ABB 

DI ,2 

; ADD 2 TO TABLE POINTER 

0020 3B01 

80 

CMP 

AX,CBX+PI] 5 COMPARTE AX REGISTER VALUE TO TABLE LOCATION 

0022 7BF9 

81 

JGE 

ILOP 

; JUMP TO INCREMENT LOOP 

0024 EB0490 

82 

00 

JMP 

INTR 

? JUMP TO INTERPOLATION ROUTINE 


bo 

84 

INTERPOLATION LOOP 


85 

DI CONTAINS ADDRESS OF THE X LOWER LIMIT VALUE 


86 

EX CONTAINS ADDRESS OF THE X UPPER LIMIT VALUE 


87 

! BP IS 

THE TEMPORARY SAVE REGISTER FOR TIE INTERPOLATION FACTOR 


88 

P 

1 



0027 33C702 

8? 

INTA: ABB 

PI, 2 

; ADD 2 TO TABLE POINTER FOR EQUfiLIZATION 


90 



: IF COMMING FROM DECREMENT LOOP 

002A 87D7 

91 

INTR: XCHG 

DX, DI 

? MOVE INDEX AND POINTER FOR WRITING OUT 

002C 8915 

92 

MOV 

CDILDX 

? UPDATE INDEX VALLE OUTSIDE LOOP 

002E 03BA 

93 

ABB 

BX.DX 

5 ADD INDEX AND BASE POINTER 

0030 8BFB 

94 

MOV 

BI,BX 

! MOVE UPPER LIMIT POINTER TO DI REGISTER 

0032 8BB0 

95 

MOV 

BX,AX 

; 'SAVE LOOK-UP VALUE IN DX REGISTER 

0034 83EF02 

96 

SUB 

DI.2 

? DECREMENT DI AND MAKE LOWER LIMIT POINTER 

0037 8B07 

97 

MOV 

AX,CBX3 

5 MOVE UPPER LIMIT TABLE VALUE TO AX 

0039 2B05 

98 

SUB 

AXtlDII 

! SUBTRACT LOWER LIMIT VALLE 

003B 8BE8 

99 

MOV 

BP, AX 

; SAVE RESULT IN BP REGISTER 

003B 2B15 

100 

SUB 

DX.CDI3 

! COMPUTE NUMERATOR OF INTERPOLATION FACTOR 

003F B80000 

101 

MOV 

AX.O 

5 CLEAR THE LEAST SIGNIFICANT BITS IN 


102 



5 DOUBLE PRECISION ACCUMULATOR TO 


303 



? PREPARE FOR DIVIDE 

0042 D1EA 

104 

SHR 

DX, 1 

1 SHIFT RIGHT 1 PLACE IN DOUBLE PRECISION 

0044 DIBS 

105 

RCR 

fiX, 1 

5 TO AVOID A DIVIDE OVERFLOW IN NEXT OPERATION 


PACE 
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MCS-P6 KACRO ASSEMBLER 

FUN 1 

3 SEPT 80 

LOG OBJ 

LINE SOURCE 

0046 D1EA 

106 

SHR DX, 1 1 SHIFT RIGHT 1 PLACE IN DOUBLE PRECISION 

0048 D1D8 

107 

RCR AX, 1 ; TO AVOIR A DIVIDE OVERFLOW IN NEXT OPERATION 

004A F7F5 

108 

DIV BP 1 DIVIDE TO COMPUTE POSITIVE 


109 

; INTERPOLATION FACTOR 

O04C 8BC8 

110 

111 ; 

MOV CX.AX 1 SAVE INTERPOLATION FACTOR IN CX REGISTER 


112 ; 

SI REGISTER INCREASED BY A FACTOR OF 2*(SI+2) 


113 ' 

! THIS ALLOWS ADDRESSING OF THE Y'S IN THE CURRENT 


114 • 

115 

116 ' 

! FUN1 TABLE OR THE NEXT FUN1A TABLE 
5 

i 

004E 8BEF 

117 

MOV BP,DI 5 SAVE LOW POINTER IN BP 

0050 83C602 

113 

ADD SI, 2 5 ADD 2 TO NUMBER OF POINTS TO COMPENSATE 


119 

5 FOR END POINTS ADDED TO DATA 

0053 D1E6 

120 

SAL SI.l 1 MULTIPLY SI BY 2 

0055 8BFE 

121 

MOV DI.SI ? SAVE NO. OF POINTS IN DI REGISTER 

0057 8800 

122 

FUN1A:M0V AX.IBX+SI] 5 MOVE UPPER Y BREAKPOINT TO AX 

0059 DIFS 

123 

SAR AX-1 ; SHIFT TO AVOID SUBTRACT OVERFLOW 


124 

; PERFORM Y INCREMENT SUBRACTIuN 

005B 3E8B12 

125 

MOV DX,BS:[BP+SU 1 MOVE LOWER Y VALUE INTO DX REG. 

005E D1FA 

126 

SAR DX, 1 ; RIGHT SHIFT 1 TO AVOID OVERFLOW IF 

127 

; BREAKPOINT SPREAD IS > 32767 

0060 2BC2 

128 

SUB AX, DX 1 COMPUTE DENOMINATOR OR INTERPOLATION FACTOR 

0062 F7E9 

129 

IMUL CX ? MULTIPLY BY X INTERPOLATION FACTOR 


130 

: SHIFT LEFT 2 TO CORRECT THE 


131 

; PREVIOUS DIVIDE SHIFT 


132 

; AND SUBTRACT SHIFTS 

0064 D1E0 

133 

SHL AX, 1 ? DOUBLE PRECISON LEFT SHIFT TO CORRECT 

0066 D1D2 

134 

RCL DX, 1 ; FOR SUBTRACT OVERFLOW SHIFTS 

0068 01 EO 

135 

136 

SHL AX,! ; DOUBLE PRECISION LEFT SHIFT TO CORRECT 

006A D1D2 

137 

RCL DX, 1 ; FOR DIVIDE OVERFLOW PREVENTION SHIFT 

006C D1E0 

138 

SHL AX, 1 ; DOUBLE PRECISION LEFT SHIFT TO CORRECT 

006E 0102 

139 

RCL DX, 1 ? FOR DIVIDE OVERFLOW PREVENTION SHIFT 

0070 3E0312 

140 

ADD DX,DS: CBP+SI] ; ADD Y LOU VALUE TO COMPLETE INTERPOLATION 

0073 03F7 

141 

ADD SI, DI 1 ADD NO. OF POINTS TO INDEX VALUE TO BE 

142 

; READY FOR FUN 1 A CALL 

0075 8BC2 

143 

MOV AX, DX 1 MOVE RESULT TO AX REGISTER FOR PL/M-86 

144 

; PROGRAM OUTPUT 

0077 C3 

145 

RET ; RETURN TO CALLING PROGRAM 

— 

146 

CODE ENDS 


147 

END 



MCS-86 MACRO ASSEMBLER FUN1 
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3 SEPT BO PAGE 4 


XREF SYMBOL TABLE LISTING 


NAME TYPE VALUE ATTRIBUTES, XREFS 

??SEG . SEGMENT SIZE=OOOOH PARA PUBLIC 

CfUN. . L NEAR OOIOH CODE 68 70# 

CfiROUP. GROUP CODE 42# 43 

CODE. . SEGMENT SIZE=0078H PARA PUBLIC -'CODE' 42# 4? 146 

HOP. . L NEAR 0014H CODE 74# 77 

FUN1. . L NEAR 0006H CODE PUBLIC 44 63# 

FUN1A . L NEAR 0057H CODE PUBLIC 44 122# 

I LOP. . L NEAR 001DH CODE 72 79# 81 

INTA. . L NEAR 0027H CODE 76 89# 

INTR. . L NEAR 002AH CODE 82 91# 

PFUN1 . L NEAR OOOOH CODE PUBLIC 44 52# 


ASSEMBLY COMPLETE, NO ERRORS FOUND 
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• SIS-! J MCS-86 HACRO ASSEMBLER V2.1 ASSEMBLY OF MODULE FUN2 
OBJECT MODULE PLACED IN :Fi:PFUN?.ORJ 

ASSEMBLER INVOKED BY: ASM86 :FJ : FEW?. SRC X'REF CATE (3 SEPT DO) 


3 SEPT go PAGE 1 


LOG OBJ 


LINE 


SOURCE 


10 

31 

12 

13 

14 

15 

16 

17 

18 
1? 
iV 
21 


*io 

*>? 


34 

^rC 

36 

•- T 
v*/ 

33 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


: FUN2 FUNCTION ROUTINE FOR A T ABLF OF TWO 

: VARIABLES WITH THE X BREAKPOINTS ALL EQUAL 

5 CA1.LNG SEQUENCE: 

? PUSH ON STACK IN THE FOLLOWING ORDER: 

: V TABLE LOCATION IN MEMORY 

: NUMBER OF X POINTS 

; X INDEX VALUE ADDRESS 

; X TABLE BASE ADDRESS 

? X VALUE TO BE LOOKED UP IN TABLE 
: PUT FOLLOWING IN REGISTERS: 

; Z INDEX VALUE ADDRESS IN PI REGISTER 
: 7 TABLF BASE ADDRESS IN BX REGISTER 

5 Z LOOK-UP VALUE IN THE AX REGISTER 

: CALL FIJN2 

; AX' REGISTER CONTAINS THE RETURNED RESULT 

i 

: PLM/86 CALL SEQUENCE 

; RSLT = FUN2(YBASE.NPTS,XINDX, XBASE, X,ZINDX,ZBASE,Z) 

? XBASE = BASE ADDRESS OF THE Y TABLE DATA 
; NPTS = NO. OF X POINTS 
; XINDX = ADDRESS OF THE X INDEX VALUE 
; XBASE = BASE ADDRESS OP THE X TABLE DATA 
; X = X LOOK-UP VALUE 
: ZINDX = ADDRESS OF THE Z INDEX VALUE 

5 ZBASE = BASE ADDRESS OF THE Z TABLE DATA 
; Z = Z LOOK-UP VALUE 

5»»»»3»»»»»»>»3 BEWARE <C<<CC<CC<<<CCC<<CCCCC<CC<<<< :«<« CCS 

T 

! THE X AND 7 TABLE INDEX LOCATIONS ARE AUTOMATICALLY UPDATED 
! BY THE ROUTINE AND THEREFORE DO NOT FOLLOW THE PL/M-86 
: BLOCK STRUCTURE 

? 

■ \ v *.*\ \*\ -• *. v % -\ -% v* \ *. \ *>. *> \ \ \ *. /- /■ r rr r r r r/ffrff r/S / .-■* ■ * ,* r r s r / ff .* ' /■ 

? j JjjtfffjJSs J Sj'j'J j/JJ.fJ s S.*J \ \SA\V. w\\ v -v. v\v v \ v\\* Cp. \.v 

7 

! ALL REGISTERS EXCEPT THE SEGMENT REGISTERS 
; ARE USED IN THE CALC. AND ARE THEREFORE DESTROYED 
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ICS-86 WORD ASSEMBLER FUN2 


3 SEPT SO PAGE 


LOG OBJ 


LINE SOURCE 


OOOO 0100 
0002 0100 
0004 0100 
0006 0100 


0010 3DFF7F 
0013 7501 

0015 48 

0016 3800 


51 

52 

53 

54 

55 

56 


59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

“ 7*7 

f / 

78 

79 
30 
81 
82 
83 
34 


NOTE: THE Y VALUES ARE ASSUMED TO 

NOT HAVE A DIFFERENCE OF GREATER THAN 32767 


FUN2A CALLNG SEQUENCES 


CALL FUN2A 

NOTE: X AND Z VALUES ARE THE SAME, Y IS DIFFERENT 

THE CX, DX AND SI REGISTERS MUST NOT BE CHANGED BEFORE THE CALL 

AX REGISTER CONTAINS THE RETURNED RESULT 


NAME FUN2 

PUBLIC FUN2,FUN2A,PFUN2 
CGROUP GROUP CODE 

DGROUP GROUP CONST, DATA, STACK, MEMORY 
ASSUME CS: CGROUP, DS: DGROUP 

DATA SEGMENT FOR THE ROUTINE RESIDES IN 
TEMPORARY STORAGE 

DATA SEGMENT PUBLIC 'BATA' 

SPC DW 1 ; STORAGE FOR PROGRAM COUNTER 

SSI DW 1 ; STORAGE FOR SI REGISTER (2 INDEX) 

YPTI DW 1 ; STORAGE FOR Y TABLE OFFSET VALUE 

DDX DW 1 ? STORAGE FOR Z INTERPOLATION FACTOR 

DATA ENDS 




85 

: BEGIN CODE SEGMENT FOR THE ROUTINE 



86 

9 



— — 


87 

CODE SEGMENT PUBLIC 

'CODE' 



88 

! UNPOP 

STACK TO 

STORE PROGRAM COUNTER 



89 

5 AND GET THE Z VARIABLES FROM THE STACK 



90 

7 



0000 3F060000 

R 

91 

PFUN2:P0P 

SPC 

? STORE PROGRAM COUNTER 

0004 58 


92 

POP 

AX 

: Z LOOK-UP VALUE 

0005 5B 


93 

POP 

BX 

! Z TABLE BASE VALUE 

0006 5F 


94 

POP 

DI 

! Z TABLE INDEX POINTER LOCATION 

0007 EB0590 


95 

JMP 

BEG 

! JUMP AROUND ASSEMBLY CALL 

OOOA 3F060000 

R 

96 

FUN2: POP 

SPC 

? SAVE PROGRAM COUNTER . 

OOOE 8B35 


97 

BEG: MOV 

SI,IDIj 

: MOVE ACTUAL INDEX VALLE INTO SI 


93 

99 

100 

101 

102 

103 

104 

105 


BEGIN C0MP6RISJ0NS TO DETERMINE Z LOCATION 


CMP AX, 32767 
JNZ NCHGZ 
DEC AX 


NCHGZ: CMP 


! CHECK IF Z VALUE IS EQUAL TO RILL SCALE 
; JUMP TO CONTINUE SEARCH LOOP 
5 DECREMENT BY ONE BIT TO PREVENT 
! Z POINTER FROM BEING LOST 
? COMPARE Z LOOK-UP VALUE WITH VALLE 


fix.cex+sn 
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KCS-86 MACRO ASSEMBLER 

FUN2 


3 SEPT 80 PAGE 3 

LOC OBJ 

LINE 

SOURCE 



106 


; IN Z TABLE POINTED TO BY LA+SI 

0013 7D0? 

107 

JGE ILOP 

; JUMP TO INCREMENTING LOOP 


103 

; DECREMENTING LOOP 


001 A 83EE02 

109 

BLOP: SUB SI, 2 

; SUBTRACT 2 FROM TABLE INDEX 

001D 3BOO 

110 

CMP AX» CBX+SI1 

: COMPARE Z LOOK-UP VALUE WITH TABLE LOS 

001F 7D0C 

111 

JGE CLZft 

! JUMP TO INTERPOLATION FACTOR COMP. 

0021 EBF7 

112 

JMP DLCP 

; JUMP TO BEGINNING OF DECREMENT LOOP 


113 

! INCREMENTING LOOP 


0023 83C602 

114 

ILOP: ADD SI, 2 

; ADD 2 TO TABLE INDEX VALUE 

0026 3B00 

115 

CMP AX,IBX+SI3 

! COMPARE Z LOOK-UP VALUE WITH TABLE LOC 

0028 7DF? 

116 

JGE ILOP 

; JUMP TO BEGINNING OF INCREMENT LOOP 

002A EB0490 

117 

JMP CLZ 

; JUHP TO INTERPOLATION FACTOR CALCULATION 


IIS 

; CLOSE Z SEARCH AND COMPUTE Z INTERPOLATION FACTOR 

002D 83C602 

119 

CLZA: ADD SI, 2 

? ADD 2 TO TABLE POINTER FOR EQUALIZATION 


120 


; IF COMMING FROM DECREMENT LOOP 

0030 03DE 

121 

CLZ: ADD BX.SI 

; ADD INDEX AND BASE Z POINTER VALLES TOGETHER 

0032 8935 

Ux. 

MOV [DI3.SI 

; SAVE Z INDEX VALLE FOR NEXT CALL 

0034 8BFB 

123 

MOV DI,BX 

; SAVE RESULT IN DI REGISTER 

0036 83EF02 

124 

SUB DI,2 

; DECREMENT POINTER FOR LOWER 7 TABLE VALUE 


125. 

; COMPUTE Z INTERPOLATION FACTOR 

0039 2B05 

126 

SUB AX,CDI3 

; SUBTRACT Z LOW FROM I INPUT VALUE 

003B 8BD0 

127 

MOV DX.AX 

; SAVE RESULT IN DX 

003D S307 

128 

KOV AX t CBX l 

; PUT Z HIGH IN AX REGISTER 

003P 2305 

129 

SUB AX.IDI3 

; SUBTRACT Z LOW FROM Z HIGH VALLE 

0041 8BC-8 

130 

MOV CX,AX 

; TEMPORARILY SAVE VALUE IN CX 

0043 3SOOOO 

131 

MOV AX,0 

: CLEAR LEAST SIG PORTION OF D I VIEWED 

0046 D1EA 

132 

SHR nx, 1 

; RIGHT SHIFT 1 TO PREVENT DIVIDE OVERFLOW 

0048 BIDS 

133 

RCR AX, 1 

: ON NEXT OPERATION 

0”j4A D1EA 

134 

SHR OX, I 

: RIGHT SHIFT 1 TO PREVENT DIVIDE OVERFLOW 

004C DIES 

135 

RCR AX,1 

; ON NEXT OPERATION 

004E F7F! 

136 

DIV CX 

; COMPUTE Z INTERPOLATION FACTOR 

0050 A30600 R 

337 

MOV CDX,AX 

; DDX NOW CONTAINS THE Z INTERPOLATION FACTOR/2 


133 

; SI, SSI CONTAIN THE 7 INDEX VALUE 

0053 89360200 R 

139 

MOV SSI, SI 

; SAVE THE Z INDEX VALUE 


140 

; POP STACK FOR VARI 

AEL.ES NECESSARY FOR X SEARCH 

0057 58 

141 

POP AX : X 

VALUE 

0058 5B* 

142 

POP BX : X 

TABLE BASF VALUE 

0059 5F 

143 

POP BI ; X 

TABLE INDEX VALUE POINTER 

005A 3835 

144 

MOV si, con : 

MOVE ACTUAL INDEX VALUE INTO SI 


145 




146 

; COMMENCE X TABLE SEARCH 

005C 3PFF7F 

147 

CMP AX, 32767 

; CHECK IF X VALUE IF EQUAL TO FULL SCALE 

005F 7501 

148 

JNZ CFUN2 

: JUMP TO CONTINUE SEARCH LOOP 

0061 48 

149 

DEC AX 

; DECREMENT BY 1 BIT TO PREVENT 


150 


! X POINTER FROM BEING LOST 

0062 3000 

151 

CFUN2: CMP AX,[BX+SD 

! COMPARE X LOOK-UP VALUE WITH VALUE 


152 


: IN X TABLE POINTED TO BY BX+SI 

0064 7009 

153 

JGE I ILOP 

; JUMP TO INCREMENTING LOOP 


154 

! DECREMENTING LOOP 


0066 83EE02 

155 

DDLOP:SUB SI, 2 

; SUBTRACT 2 FROM TABLE INDEX 

0069 '3F00 

156 

OP AX.CBX+SII 

; COMPARE X LOOK-UP VALUE WITH TABLE LOC 

006B 7D0C 

157 

JGE CMPA 

; .JUMP TO INTERPOLATION FACTOR COMP. 

006D EPF7 

15S 

.JMP DDLOP 

: JUMP TO BEGINNING OF DECREMENTING LOOP 


159 

; INCREMNTING LOOP 


006F 33C602 

160 

I ILOP: ADD SI, 2 

; ADD 2 TO TABLE INDEX VALUE 
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MCS-36 MACRO ASSEMBLER 

FUN2 



3 SEPT 

LOG OBJ 

LINE SOURCE 



0072 3E00 

161. 

CMP 

AX.PBX+SI3 

COMPARE X LOOK-UP VALUE WITH TABLE LOC. 

0074 7DF9 

162 

JGE 

IILOP 

: JUMP TO INCREMENTING LOOP 

0076 EB0490 

163 

..IMP 

CMPP 

JUMP TO INTERPOLATION FACTOR CALCULATION 


164 ! 

: CLOSE X 

SEARCH AND PREPARE TO COMPUTE X INTRP FACTOR 

0079 S3C602 

165 CMP A: ADD 

nr n « 

v*l Ji 1 

ADD 2 TO TABLE POINTER FOR EQUALIZATION 


166 


1 

IF COMMING FROM DECREMENT LOOP 

007C 03DE 

167 CHRP: ADD 

BX.SI 1 

ADD INDEX AND TABLE POINTER BASE VALUES 

007E 3935 

163 

MOV 

CDII.SI i 

SAVE INDEX POINTER VALLE IN SPEC. LOC. 

0030 8BFB 

169 

MOV 

DIiBX ! 

SAVE RESULT IN DI REGISTER 

0032 33EF02 

170 

SUB 

DIi 2 1 

DECREMENT POINTER FOR LOWER X TABLE VALUE 


1.71 ! 

: COMPUTE X INTERPOLATION FACTOR 

0035 2B05 

172 

SUB 

AX, [DII 

SUBTRACT X LOW FROM X INPUT VALUE 

0087 8BD0 

173 

MOV 

DX,AX ! 

SAVE RESULT IN DI 

0039 8B07 

174 

MOV 

AXiCBXI 

PUT X HIGH IN AX REGISTER 

0088 2805 

175 

SUB 

AX, [DII ! 

SUBTRACT Z LOW FROM Z INPUT VALUE 

008D 3BC8 

176 

MOV 

CX,AX 

SAVE VALUE IN CX 

008F 880000 

177 

MOV 

AX,0 { 

CLEAR LEAST SIGNIFICANT PORTION OF DIVIDEND 

0092 D1EA 

178 

SHR 

DX , 1 

SHIFT RIGHT 1 PLACE IN DOUBLE PRECISON 

0094 D1D8 

179 

RCR 

AX, 1 ' 

TO AVOID A DIVIDE OVERFLOW IN NEXT OPERATION 

0096 D1EA 

130 

SHR 

DX, 1 

SHIFT RIGHT 1 PLACE IN DOUBLE PRECISION 

0093 D1D8 

181 

RCR 

AX,1 i 

TO AVOID A DIVIDE OVERFLOW IN NEXT OPERATION 

009A F7F! 

182 

DIV 

CX 

DIVIDE TO COMPUTE INTERPOLATION FACTOR 

009C SBC8 

183 

MOV 

CX,AX i 

SAVE RESULT IN CX REGISTER 


184 • 

f 




185 ' 

; CX CONTAINS THE X INTERPOLATION FACTOR/2 


186 ! 

1 DDX CONTAINS THE Z INTERPOLATION FACTOR/2 


137 i 

» 




188 ! 

: SAVE THE INDEX REGISTER VALUES SO 


189 ! 

! SO THAT THEY CAN BE RECALLED LATER FOR 


190 ! 

1 PERMANENT STORAGE OUTSIDE ROUTINE 


191 

! 




192 ! 





193 ! 

1 COMPUTE THE Y TABLE OFFSETS FROM THE 


194 i 

; START OF THE Y TABLE 


009E 58 

195 

POP 

AX 

AX CONTAINS THE NUMBER OF X F'TS. 

009F 050200 

196 

ADD 

AX, 2 

ADD 2 TO NO. OF X POINTS FOR END POINTS 

00A2 8BE8 

197 

MOV 

BP, AX 

SAVE RESULT IN BP 

00A4 DIES 

193 

SHL 

BP, 1 

MULTIPLY RESULT BY 2 

OOA6 F7260200 R 

199 

MUL 

SSI ! 

MULTIPLY BY Z INDEX VALLE 

OOAA 03C6 

200 

ADD 

AX, SI 

ADD X INDEX VALUE 




201 

' 




202 

• 

? 

AX CONTAINS THE OFFSET OF THE HIGH 



203 

> 

Y WORD FROM THE START OF THE Y TABLE 



204 

205 
2X56 

• 

y 

OFFSET (BYTES) = Z OFFSET * (X POINTS +2) + X OFFSET 




STORE TOTAL NUMBER OF X POINTS + 2 IN BYTES 

OOAC 8BF5 


207 


MOV SI, BP 



208 

• 

STORE Y TABLE OFFSET FOR POSSIBLE FUN2A CALL 

OOAE A30400 

R 

209 


MOV YPTI.AX 

GOBI EB0390 


210 


JMP CFUN ; JUMP TO CONTINUE FUNN2 ROUTINE 



211 

• 




212 

• 

ENTRY POINT FOR THE FUN2A CALL 



213 

■ 


0084 53 


214 

FUN2ASP0P AX 

GOES A30000 

R 

215 


MOV SPC, AX ; POP AND STORE PROGRAM COUNTER 


PAGE 4 


5 POP AND STORE PROGRAM COUNTER 



22 


MCS-86 MACRO ASSEMBLER 


3 SEPT 80 PAGE 


LOG 

OBJ 


LINE 

SOURCE 


GOBS 

A) 0400 

R 

216 

►—i 

h— 

Ct_ 

X 

<L 

z> 

5 RESTORE Y TABLE OFFSET TO AX 




217 






218 

; COMPUTE BASE POINTER VALUES USING Y LOACTION 


219 ; AND THE COMPUTED OFFSET 

220 ; 


OOBB CD 

221 

CFUN: POP 

BP 

5 BP CONTAINS Y TABLE BASE VALUE 

ODBC 8RF3 

222 

MOV 

DIi AX 

: MOVE Y TABLE OFFSET TO DI . 

OOBE 03FD 

223 

ADD 

DI » BP 

5 ADD Y TABLE BASE VALUE 

wCu 8BD- 

224 

MOV 

BX.BI 

; MOVE RESULT TO BX REGISTER 

00C2 83EB0? 

225 

SUB 

BX,2 

; SUBTRACT 2 FROM RESULT 


2x.6 

/ 

1 

i DI AND 

BX' CONTAIN THE HIGH AND LOW TABLE 



! LOCATION VALUES FOR 

INTERPOLATION 



; SI CONTAINS NUMBER 

OF X POINTS 


230 

1 




231 

! COMPUTE THE Y HIGH 

INTERPOLATED VALUE 

OOC5 8BQ5 

232 

MOV 

AX.CDI1 

; MOVE Y HIGH PRIMITIVE VALUE 

OOC7 DIES 

233 

SAR 

AX.l 

j SHIFT VALUE TO AVOID OVERFLOW 


234 



; IF BREAKPOINT SPREAD >32768 

OOC'7 SP2F 


MOV 

EP.CBX] 

; MOVE Y LOW FRIMITVE VALUE 

OOCB D1FD 

236 

SAR 

BP, 1 

! SHIFT VALUE TO AVOID OVERFLOW 

OOCD 2BC5 

x-0/ 

SUB 

AX, BP 

? SUBTRACT Y LOW PRIMITIVE VALUE 

OOCF F7E9 

OO O 

IMUL 

CT 

? MULTIPLY X INTERPOLATION FACTOR 

00D1 D1E0 

onn 

Av'“ 

SHL 

AX, 1 

; SHIFT RESULT TO CORRECT FOR 

00D3 B1B2 

240 

RCL 

DX,1 

5 SUBTRACTION OVERFLOW SHIFT 

00D5 D1E0 

241 

'SHL 

ftV * 

HA> i 

; MULTIPLY RESULT BY 2 TO CORRECT FOR 

0007 D1D2 

242 

RCL 

DX.l 

; A PREVIOUS DIVIDE SHIFT 

non? DIED 

243 

SHL 

AX.l 

; MULTIPLY RESULT BY 2 TO CORRECT FOR 

OODB D1D2 

244 

RCL 

DX , 1 

; A PREVIOUS DIVIDE SHIFT 

OODD 0317 

245 

ADD 

DX.tBXl 

; ADD Y LOW PRIMITIVE VALUE 

OOBF 8BEA 

246 

MOV 

BP,DX 

; BP NOW CONTAINS THE Y HIGH INTERPOLATED VAL! 


iAi 

24.-: 

; COMPUTE THE NEW INDEX VALUES FOR Y LOW 


249 

; INTERPOLATED VALUE COMPUTATION 

OOE1 2BFE 

250 

SUB 

DI .SI 

! COMPUTE NEW Y HIGH PRIMITIVE VALUE LOCATION 

00E3 2BDE 

251 

SUB 

BX, 51 

; COMPUTE NEW Y LOW PRIMITIVE VALUE LOCATION 


252 

? COMPUTE THE Y LOW 

INTERPOLATED VALUE 

00E5 8B05 

253 

MOV 

AX, EDI3 

; MOVE Y HIGH PRIMITIVE VALUE 

00E7 DIFS 

254 

SAR 

AX.l 

; SHIFT VALUE TO AVOID OVERFLOW 


-“.ETC* 

J-J-J 



? IF BREABPOINT SPREAD I > 32767 

00E9 8B37 

256 

MOV 

SI.CBX3 

; MOVE Y LOW PRIMITIVE VALUE 

OOEB D1FE 

257 

SAR 

SI,1 

; SHIFT VALUE TO AVOID OVERFLOW 

OOED 2BC6 

258 

SUB 

AX, SI 

; SUBTRACT Y LOW PRIMITIVE VALUE 

OOEF F7E9 

259 

IMUL 

CX 

: MULTIPLY X INTERPOLATION FACTOR 

OOF1 DIED 

260 

SHL 

AX.l 

5 SHIFT RESULT TO CORRECT FOR 

OOFS D1D2 

261 

RCL 

DX, 1 

5 SUBTRACTION OVERFLOW SHIFT 

OOFS D1E0 

262 

SHL 

AX.l 

; MULTIPLY RESULT BY 2 TO CORRECT FOR 

OOF7 D102 

263 

RCL 

DX.l 

; A PREVIOUS DIVIDE SHIFT 

OOFS D1E0 

264 

SHL 

AX.l 

; MULTIPLY RESULT BY 2 TO CORRECT FOR 

OOFB D1D2 

265 

RCL 

DX.l 

5 A PREVIOUS DIVIDE SHIFT 

OOFD 0317 

266 

ADD 

DX, [BX] 

? ADD Y LOW PRIMITIVE VALUE 

OOFF SBF2 

267 

MOV 

SI.DX 

; SI CONTAINS Y LOW INTERPOLATED VALUE 


268 

? 




269 

; COMPUTE THE FINAL 

INTERPOLATED V VALUE 

0101 83C5 

270 

MOV 

AX, BP 

; MOVE Y HIGH INTERPOLATED VALLE 



23 



3 SEPT 80 PAGE 6 


LOG OBJ LINE SOURCE 


0103 DIFS. 

271 

SAR 

AX.l ; 

SHIFT TO AVOID SUBTRACT OVERFLOW 

0105 PIFfl 

272 

SAR 

DX,1 ! 

SHIFT TO AVOID SUBTRACT OVERFLOW 

0107 2BC2 

273 

SUB 

AX.DX i 

SUBTRACT Y I.OW INTERPOLATED VALUE 

010? F72E0600 

R 274 

IMUL 

DDX ! 

MULTIPLY BY Z INTERPOLATION FACTOR 

010D D1E0 

275 

SHL 

AX.l ! 

SHIFT RESULT TO CORRECT FOR 

01 OF D1D2 

276 

RCL 

DX.l ' 

SUBTRACTION OVERFLOW SHIFT 

0111 D1E0 

277 

SHL 

AX.l ' 

MULTIPLY BY 2 TO CORRECT FOR 

0113 D1D2 

278 

RCL 

DX.l ' 

A PREVIOUS DIVIDE SHIFT 

0115 D1E0 

27? 

SHL 

AX.l i 

MULTIPLY BY 2 TO CORRECT FOR 

0117 D1D2 

280 

RCL 

DX.l 

A PREVIOUS DIVIDE SHIFT 

0119 03D6 

231 

ADD 

DX.SI • 

ADD Y LOW INTERPOLATED VALUE 


282 ? PREPARE TO RETURN TO THE MAIN PROGRAM 

283 ; PUT PROGRAM COUNTER BACK ON STACK 

01 IB 3B2EOOOO R 284 MOV BP.SPC 

01 IF 55 285 PUSH BP 

286 i PUT RESULT IN AX REGISTER FOR PROPER TRANSFER 

0120 8BC2 237 MOV AX.BX 

0122 C3 283 RET 

23? CODE ENDS: 

290 END 
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3 SEPT SO PAGE 7 


MCS-86 MACRO ASSEMBLER FUN2 


VREF SYMBOL TABLE LISTING 


NAME 

TYPE 

VALUE 

ATTRIBUTES, XREFS 


?''SEG . 

SEGMENT 


SIZE=OOOOH PARA PUBLIC 


BEG . . 

L NEAR 

OOOEH 

COPE 95 97# 


CFUN. . 

L NEAR 

OOBBH 

CODE 210 221# 


CFUN2 . 

L NEAR 

0062H 

CODE 14S 151# 


CGROUP, 

GROUP 


COPE 71# 73 


CLZ . . 

L NEAR 

0030H 

CODE 117 121# 


CL7A. . 

L NEAR 

002DH 

CODE 111 119# 


cm . 

L NEAR 

0079H 

CODE 157 165# 


CMPP. . 

1. NEAR 

007CH 

CODE 163 167# 


COPE. . 

SEGMENT 


SIZE=0123H PARA PUBLIC ' 

CODE' 71# 37 28' 

CONST . 

SEGMENT 


SIZE=OOOOH —UNDEFINED— 

72# 

DATA. . 

SEGMENT 


SIZE=0008H PARA PUBLIC ' 

DATA' 72# 78 33 

DDL OP . 

L NEAR 

0066H 

CODE 155# 15S 


DDK . . 

V WORD 

0006H 

DATA 82# 137 274 


DGROIJP. 

GROUP 


CONST DATA STACK MEMORY 

72# 73 

PLOP. . 

L NEAR' 

001 AH 

CODE 109# 112 


FIJN'2. . 

L NEAR 

OOOAH 

CODE PUBLIC 70 96# 


FUN2A . 

L NEAR 

00B4H 

CODE PUBLIC: 70 214# 


I I LOP . 

L NEAR 

006FH 

CODE 153 160# 162 


ILOP. . 

L NEAR 

0023H 

CODE 107 114# 116 


MEMORY. 

SEGMENT 


SIZE=OOOOH -UNDEFINED- 


NCHGZ . 

L NEAR 

0016H 

CODE 102 105# 


PFUN2 . 

L NEAR 

OOOOH 

CODE PUBLIC 70 91# 


SPG . . 

V WORD 

OOOOH 

DATA 79# 91 % 215 234 


SSI . . 

V WORD 

0002H 

DATA 80# 139 199 


STACK . 

SEGMENT 


8!ZE=0000H -UNDEFINED- 


YPTI. . 

V WORD 

OOOAH 

DATA 31# 209 216 



ASSEMBLY COMPLETE, NO ERRORS FOUND 
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HCS-36 MACRO ASSEMBLER FUNS 


3 SEPT SO PAGE 1 


ISIS-II MCS-86 MACRO ASSEMBLER V2.1 ASSEMBLY OF MODULE FUNS 
OBJECT MODULE PLACED IN :F1 :PFUN3.0BJ 

ASSEMBLER INVOKED BY: ASMS*. : FI : PFUN3. SRC XREF DATE (3 SEPT 80) 


LOG OBJ 


LINE SOURCE 


1 


4 

5 

6 

7 
S 

8 
10 
11 
12 

13 

14 

15 

16 
17 
13 

19 

20 
21 
22 

23 

24 

'75 

26 

27 

28 
•*v% 

l y 

30 

31 

32 
00 

34 

35 

36 

37 

33 

39 

40 
43 

42 

43 

44 

45 

46 

47 

48 

49 

50 


FUNS FUNCTION ROUTINE FOR TWO VARIABLE TABLES 
X BREAKPOINTS NOT EQUAL 


CALLNG SEQUENCE: 


PUSH ON THE STACK IN THE FOLLOWING ORDER 
Y TABLE BASE VALLE 
X LOOK-UP VALUE 
X TABLE BASE ADDRESS 
NUMBER OF X POINTS 
PUT FOLLOWING IN REGISTERS 
Z TABLE INDEX VALUE ADDRESS IN DI REGISTER 
Z TABLE BASE ADDRESS IN BX REGISTER 
Z LOOK-UP VALUE IN AX REGISTER 
CALL FUN3 


AX REGISTER CONTAINS THE RETURNED RESULT 
Z VALLE MUST BE LIMITED TO THE Z TABLE RANGE 




PL/M-36 CALLING SEQUENCE 


RSLT = FUN3(YBASE,X, XBASE, NPTS,ZINDX,ZBASE,Z) 


YBASE = V TABLE EASE ADDRESS 

X = X LOOK-UP VALUE 

XBASE = BASE ADDRESS OF THE X TABLE 

NOPTS = NO. OF X POINTS IN TABLE 

ZINDX = ADDRESS OF THE Z INDEX VALUE 

ZBASE = BASE ADDRESS OF THE Z TABLE 

Z =.Z LOOK-UP VALUE 


»»7}»}'>VX>>YV» 


,\\\\\\ 
j .• j j j / 


BEWARE <<CC<<<CC<C<C<<<<<<<<<<<<<<<<<<<<<< 


THE Z TABLE INDEX LOCATION IS AUTOMATICALLY UPDATED 
' BY THE ROUTINE AND THEREFORE DO NOT FOLLOW THE PL/M-86 
BLOCK STRUCTURE RULES 

>»»»»»»»»»»»»»»»C«C««««««««C««««<««CC<« 


fiiHmimHmHmHHHllllHHiHlHIHmtHHHHiHtlltHHIl 
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MCS-36 MACRO ASSEMBLER 

FUN3 



3 SEPT & 

LOC OBJ 

LINE 

SOURCE 




51 

7 




52 

NAME FUN3 




53 

PUBLIC FUN3.PFUN3 




54 

CGRCUP GROUP CODE 




55 

DOROUP GROUP DATA, CONST, STACK, MEMORY 


56 

ASSUME CS: CGROUP, DS: DGROUP 




57 

7 




58 

; DATA SEGMENT 



— 

59 

DATA SEGMENT PUBLIC 'DATA'' 



OOOO 0100 

60 

SPC3 DM 3 

: 

STORE PROGRAM COUNTER 

0002 0100 

61 

DDXZ DM 1 

7 

Z INTERPOLATION FACTOR 

0004 0100 

62 

NBYT DU 1 

7 

NUMBER OF X BYTES PER Z VALLE 

0006 0100 

63 

XSAV DU 1 

7 

X LOOK-UP VALUE 

0008 0100 

64 

HXOFF DU 1 

« 

7 

HIGH X OFFSET 

OOOA 0100 

65 

LXOFF DM 1 

* 

7 

LOM X OFFSET 

OOOC 0100 

66 

DDX DU 1 

j 

X INTERPOLATION FACTOR 

OOOE 0100 

67 

XB DM 1 

» 

7 

BREAKPOINT COORDINATE 

0010 0100 

68 

XC DM 1 

■ 

BREAKPOINT COORDINATE 

0012 0100 

6? 

XE DM 1 

« 

7 

BREAKPOINT COORDINATE 

0014 0100 

70 

YB DM 1 

; 

BREAKPOINT COORDINATE 

0016 0100 

71 

YC DU 1 

; 

BREAKPOINT COORDINATE 

0018 0100 

)i. 

YE DM 1 

; 

BREAKPOINT COORDINATE 

— 

. 73 

Tfl 

DATA ENDS 

* 




/4 

75 

? CODE SEGMENT 



— 

76 

CODE SEGMENT PUBLIC "CODE' 




77 

7 




78 

! START FUNCTION ROUTINE 




79 

• 



0000 8F060000 R 

80 

PFUN3: POP SPC3 


STORE PROGRAM COUNTER 

0004 58 

31 

POP AX 

7 

Z LOOK-UP VALLE 

0005 5B 

on 

POP BX 

7 

Z TABLE BASE VALUE 

0006 5F 

83 

POP BI 

7 

Z TABLE INDEX VALUE ADDRESS 

0007 EB0590 

84 

JMP BEG 

7 

JUMP ASSEMBLY LANGUAGE CALL 

OOOA 8F060000 R 

85 

FUNS ! POP SPC3 

1 

SAVE PC FOR ASSEMBLY CALL 

OOOE 8B35 

86 

BEG: MOV SI, ton 

• 

« 

LOAD Z TABLE INDEX VALUE POINTER 


37 

1 




38 

oo 

; DETERMINE Z LOCATION 

• 



0010 3B00 

07 

90 

CMP AX.CBX+SI3 

• 

COMPARE I LOOK-UP VALUE MITH VALLE IN 


91 


7 

Z TABLE POINTED TO BY BX+SI 

0012 7D09 

92 

JOE IL0P3 

7 

JUMP TO INCREMENTING LOOP 


93 





94 

! DECREMENTING LOOP 




95 

• 



0014 83EE02 

96 

DL0P3: SUB SI, 2 

■i 

SUBTRACT FROM TABLE INDEX 

0017 3B00 

97 

CMP AX,[BX+SI3 

• 

7 

COMPARE Z LOOK-UP VALUE MITH TABLE LOC. 

001? 7D0C 

98 

JGE CLZA3 

• 

7 

JUMP TO INTERPOLATION FACTOR CALC 

00 IB EBF7 

99 

JMP DL0P3 

> 

JUMP TO BEGINNING OF DECREMENT LOOP 


100 

7 




101 

! INCREMENTING LOOP 




102 

7 



00 ID 83C602 

103 

IL0P3: ADD SI, 2 


ADD 2 TO THE INDEX POINTER 

0020 3B00 

104 

CMP AX, CBX+SI3 

• 

1 

COMPARE Z LOOK-UP VALUE MITH TABLE LOC. 

0022 7DF9 

105 

JGE IL0P3 

; 

JUMP TO BEGINNING OF INCREMENT LOOP 


PAGE 2 


MCS-86 MACRO ASSEMBLER 

FUNS 



I.OC OBJ 

LINE 

SOURCE 


0024 EB04?0 

106 

•JMP 

CLZ3 


107 

; 



103 

; CLOSE z 

SEARCH AND 


10? 

7 


0027 83C602 

110 

CLZA3: ADD 

SI, 2 


113 



002A 03DE 

112 

CLZ3: ADD 

BX.SI 

002C 8935 

113 

MOV 

IDIL SI 

002E 8BFB 

114 

MOV 

DIiBX 

0030 83EF02 

115 

SUB 

DI,2 

0033 2B05 

116 

SJJB 

AX, [DII 

0035 8BD0 

117 

MOV 

DX , AX 

0037 8B07 

113 

MOV 

AXiCBX] 

0039 2B05 

239 

SUB 

AXiIDII 

003B sees 

120 

MOV 

CXiftX 

003D BSOOOO 

121 

MOV 

AX,0 

0040 D1EA 

122 

SHR 

DX,3 

0042 D1DS 

123 

RCR 

AX. 1 

0044 D1EA 

124 

SHR 

DX, 1 

0046 D1D8 

125 

RCR 

AX,1 

0043 F7F! 

126 

DIV 

CX 

004A A30200 R 

127 

Moy 

DDXZ, AX 


i nn 




liO 

? 



12? 

: COMPUTE OFFSET AND A 


130 

: X OFFSET = (4 OF 


131 

1 NBYT = 

(# OF XPOI 

004D 53 

332 

FOP 

AX 

004E 050200 

133 

ADD 

AX, 2 

0051 8BE3 

334 

MOV 

BP, AX 

0053 D1E5 

135 

SHL 

BP, 3 

0055 F7E6 

136 

MUL 

SI 

0057 8BF0 

137 

MOV 

SI, AX 

0059 892E0400 R 

138 

MOV 

NBYT, BP 


IS 1 ? 




140 

; 



141 

? COMPUTE X 

LOCATION F 

005D 5F 

142 

POP 

DI 

005E 8BD8 

143 

MOV 

BX, AX 

0060 53 

144 

POP 

AX 
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3 SEPT 80 


JUMP TO INTERPOLATION LOOP 


SEARCH AND COMPUTE Z INTERPOLATION FACTOR 


145 

146 

147 

148 
14 ? 


ADD 2 TO TABLE INDES FOR EQUALIZATION 

IF COMING FROM DECREMENT LOOP 

ADD INDEX AND BASE Z POINTER VALUES 

SAVE Z INDEX VALUE IN PROPER ADDRESS 

SAVE RESULT IN 01 REGISTER 

DECREMENT POINTER FOR LOWER Z TABLE VALUE 

SUBTRACT Z LOW FROM Z INPUT VALUE 

SAVE RESULT IN DX 

PUT Z HIGH IN AX REGISTER 

SUBTRACT Z LOW FROM Z HIGH VALUE 

TEMPORARILY SAVED VALUE IN CX 

CLEAR LEAST SIG. PORTION OF DIVIDEND 

FULL RIGHT SHIFT TO PREVENT OVERFLOW 

ON THE NEXT DIVIDE 

FULL RIGHT SHIFT TO PREVENT OVERFLOW 

ON THE NEXT DIVIDE 

DIVIDE TO COMPUTE Z INTERPOLATION FACTOR 
DDXZ CONTAINS THE Z INTERPOLATION FACTOR 


AX CONTAINS THE NUMBER OF X POINTS 
ADD 2 TO ACCOUNT FOR END POINTS 
SHIFT TO BP REG. 

MUL. BY 2 TO GET BYTES/Z VALUE 
MIL BY Z INDEX TO GET TOTAL OFFSET 
SI CONTAINS THE HIGH X OFFSET 
NBYT CONTAINS THE NUMBER OF X BYTES 
PER Z VALLE (VECTOR) 


DI CONTAINS X TABLE BASE VALUE 
EX CONTAINS THE HIGH X OFFSET 
X LOOK-UP VALUE 


AX CONTAINS THE X LOOK-UP VALUE 
BP CONTAINS THE NUMBER OF X BYTES PER VECTOR 
DI CONTAINS THE X TABLE BASE VALUE 


0061 3DFF7F 


150 

151 

: DETERMINE X 
CMP 

OFFSET FROM 
AX, 32767 

0064 7501 


152 

JNZ 

CFUN3 

0066 43 


153 

DEC 

AX 

0067 A30600 

R 

154 

155 

7 

CFUN3: MOV 

XSAV,AX 

006A 3B01 


156 

CMP 

AX, CDI+BX] 

006C 7D0? 


157 

158 

JGE 

ILPXH3 


35 ? 

160 


CHECK IF X VALUE IS PLUS FULL SCALE 
JUMP TO CONTINUE SEARCH ROUTINE 
DECREMENT X VALUE TO PREVENT POINTER 
FROM GETTING LOST 
SAVE X LOOK-UP VALUE 
COMPARE X LOOK-UP VALUE WITH VALUE 
IN TABLE POINTED TO BY DI+BX 
JUMP TO INCREMENTING LOOP 


DECREMENTING LOOP 


PAGE 



MCS-86 MACRO A? 

SEMBLER 

FUNS 


28 

3 SEPT 80 PAGE 4 

LOC OBJ 


LINE 

SOURCE 



006E 83EB02 


161 

DLPXH3:SUB 

BXi2 

: SUBTRACT 2 FROM INDEX VALUE 

0071 3B01 


162 

CMP 

AX.CDI+BX3 

5 COMPARE X LOOK-UP VALUE TO TABLE 

0073 7D0C 


163 

JOE 

CLXA3 

5 JUMP TO INTERPOLATION LOOP 

0075 EBF7 


164 

JMP 

DLPXH3 

; JUMP TO BEGINNING OF DECREMENT LOOP 



165 

i 

? 





166 

! INCREMENTS 

LOOP 


0077 S3C302 


167 

ILPXH3: ADD 

BX.2 

; ADD 2 TO INDEX POINTER 

007A 3B01 


168 

CMP 

AX,[DI+RX] 

: COMPARE X LOOK-UP VALUE TO TABLE 

007C 7DF9 


169 

JOE 

ILPXH3 

; JUMP TO INTERPOLATION LOOP 

007E EB0490 


170 

JMP 

CLX3 

; JUMP TO BEGINNING OF INCREMENTING LOOP' 



171 

0 





172 

! FINAL COMPUTATION OF FINAL X LOCATIONS FOR Z HIGH VALUE 



173 

7 



0081 83C302 


174 

CLXA3: ADD 

BX.2 

INCREMENT POINTER BY 2 IF 



175 



COMING FROM DECREMENT LOOP 

0084 8BEF 


176 

CLX3: MOV 

BP.DI 

BP CONTAINS X TABLE BASE VALUE 

0086 891E0800 

R 

177 

MOV 

HXOFF.BX 

HXOFF CONTAINS THE X HIGH OFFSET (XE) 

008A 03FB 


178 

ADD 

01, BX 


003C 8BDF 


179 

MOV 

BX.DI 

DI HAS LOCATION OF HIGH X (XE) FOR Z HIGH 

008E 83EB02 


180 

SUB 

BX.2 

BX HAS LOCATION OF LOW X (XB) FOR Z HIGH 

0091 893E1200 

R 

181 

MOV 

XE.DI 

SAVE LOCATION IN XE 

0095 891E0E00 

R 

182 

183 

134 

MOV 

XB.BX 

SAVE LOCATION IN XB 



7 

5 COMPUTE LOW X LOCATIONS FOR LOW VALUE OF Z 



135 

* 





186 

; COMPUTE LOW X LOCATION 




187 

; X OFFSET 

= (X OFFSET - NBYT) 

0099 2B360400 

R 

188 

SUB 

SI, NBYT : 

SI CONTAINS LOW X VECTOR OFFSET 

009D 8BDE 


139 

MOV 

BX.SI 

BX CONTAINS X VECTOR OFFSET FOR Z LOW 

009F 8BFD 


190 

MOV 

DI.BP 

SI CONTAINS X TABLE BASE VALUE 



191 

; DETERMINE X LOCATION 


OOA1 A 10600 

R 

192 

MOV 

AX.XSAV : 

AX CONTAINS X LOOK-UP VALUE 

OOA4 3B01 


193 

CMP 

AX.CDI+BX] 

COMPARE X LOOK-UP VALLE WITH VALUE 



194 



IN X TABLE POINTED TO BY DI+BX 

00A6 7B09 


195 

JOE 

ILPXL3 

JUMP TO INCREMENTING LOOP 



196 

; DECREMENTING LOOP 


00A3 83EB02 


197 

DLPXL3JSUB 

BX,2 

DECREMENT INDEX POINTER 

OOAB 3B01 


193 

CMP 

AX, CDI+BX3 

COMPARE X LOOK-UP VALUE WITH TABLE LOC. 

OOAD 7B0C 


199 

JGE 

CLXLA3 

JUMP TO FINAL X LOW CALC. 

OOAF EBF7 


200 

JMP 

DLPXL3 

JUMP TO DECREMENT LOOP 



201 

7 





202 

5 INCREMENTING LOOP 


00B1 83C302 


203 

ILPXL3: ADD 

BX , 2 

INCREMENT INDEX POINTER 

00B4 3B01 


204 

CMP 

AX.TDI+BXJ 

COMPARE X LOOK-UP VALUE WITH TABLE LOC. 

00B6 7BF9 


205 

JGE 

ILPXL3 

JUMP TO INCREMENT LOOP 

OOB8 EB0490 


206 

JMP 

CLXL3 

JUMP TO FINAL X LOW CALC. 



207 

7 





203 

i COMPUTE X LOW LOCATIONS FOR LOW VALUE OF Z 



209 

; X OFFSET = 

(X OFFSET-NBYT) 


OOBB 83C302 


210 

CLXLA3J ADD 

BX,2 

ABB 2 TO INDEX POINTER FOR EQUALIZATION 



211 



IF COMMING FROM DECREMENT LOOP 

OOBE 03FB 


212 

CLXL3: ADD 

DI.BX 

ADD OFFSET AND BASE POINTERS 

OOCO 391E0A00 

R 

213 

MOV 

LXOFF, BX 

LXOFF CONTAINS THE LOW X OFFSET (XD) 

00C4 3BDF 


214 

MOV 

BX.DI 

DI CONTAINS LOCATION OF HIGH X (XD) FOR 7 LOW 

00C6 83EB02 


215 

SUB 

BX,2 

BX CONTAINS LOCATION OF LOW X (XC) FOR Z LOW 


29 


HCS-36 MACRO A 

SSEMBLER 

FUNS 




3 SEPT 80 

LOC OBJ 


LINE 

SOURCE 




OOP? 891E1000 

R 

216 

217 

218 

219 

220 

MOV 

XC.BX 

" 

SAVE LOCATICN IN XC 



; COMPUTE XG 

= XD+DDXZ1XE-XD) 



OOCD 8B1E1200 

R 

MOV 

EX' » XE 

’ 

BX CONTAINS HIGH X LOCATION 

OOD1 8B07 


221 

MOV 

AX.CBXI 

7 

MOVE XE TO AX 

00D3 D1F8 


222 

SAR 

AX. 1 

• 

7 

RIGHT SHIFT 1 IN CASE POINT 



O'W 



* 

7 

DELTA IS GREATER THAN 32767 

00D5 8B15 


224 

MOV 

px,crn 

* 

■l 

MOVE XD 

OOB7 D1FA 



SAR 

DXv 1 

7 

SHIFT TO AVOID SUBTRACT OVERFLOW 

00D9 2BC2 


226 

SUB 

AX,DX 

* 

7 

SUBTRACT XD 

GODS F72E0200 

R 

227 

IHUL 

DDXZ 

« 

7 

MULTIPLY BY INTERPOLATION FACTOR 

OODF D1EO 


223 

SHL 

AXvl 

7 

SHIFT TO CORRECT FOR 

OOE1 D1D2 


229 

RCL 

DX» 1 

7 

SUBTRACT OVERFLOW SHIFT 

00E3-D1E0 


230 

SHL 

AXvl 

* 

7 

MULTIPLY RESULT BY 2 

00E5 D1B2 


'231 

RCL 

DXvl 

• 

7 

TO MAKE UP FOR SHIFT BEFORE DIVIDE 

OOE7 OlEO 


232 

SHL 

AX,1 

7 

FULL LEFT SHIFT TO CORRECT 

00E9 B1D2 


233 

RCL 

DX,1 

7 

PREVIOUS SHIFT 

OOEB 0315 


n~t a 
jLO* t 

ADD 

DX,CDI3 

* 

1 

ADD XD TO RESULT 

OOED 3BEA 


c 

■fc-WJ-J 

! 

MOV 

BP, DX 

7 

BP CONTAINS XG 



237 

; COMPUTE XF 

=XC+DDXZ ( XB-XC) 





238 





OOEF 8B1E0E00 

R 

wc:o7 

MOV 

PX , XB 


BX CONTAINS X HIGH LOCATION 

00F3 SEE 1000 

R 

240 

MOV 

DI,XC 

7 

DI CONTAINS X LOW LOCATION 

00F7 8B07 


241 

MOV 

AX.IBX] 


MOVE XE TO AX REGISTER 

00F9 D1F3 


242 

SAR 

AX, 1 

7 

RIGHT SHIFT 1 IN CASE POINT 



243 




SPREAD IS GREATER THAN 32767 

OOFB 8B15 


244 

MOV 

DX, CPI] 

7 

MOVE XC TO DX REG. 

OOFD D! FA 


245 

SAR 

DX, 1 


SHIFT TO AVOID SUBTRACT OVERFLOW 

OOFF 2BC2 


246 

SUB 

AX, BX 

> 

SUBTRACT XC 

0101 F72E0200 

R 

247 

IHUL 

DDXZ 

7 

MULTIPLY BY INTERPOLATION FACTOR 

0105 D1E0 


248 

SHL 

AX,1 

7 

SHIFT TO CORRECT FOR 

0107 D1D2 


249 

RCL 

DX, 1 


SUBTRACT OVERFLOW SHIFTS 

0109 D1E0 


250 

SHL 

AX,1 

7 

MULTIPLY BY 2 TO MAKE UP FOR 

010B D1D2 


251 

RCL 

DX, 1 


PREVIOUS SHIFT BEFORE DIVIDE 

010D B1E0 


252 

SUL 

AX,1 

7 

FULL LEFT SHIFT TO CORRECT 

010F D1D2 


253 

RCL 

DX,! 

• 

j 

PREVIOUS SHIFT 

0111 0315 


254 

ADD 

DX.CDI] 

■ 

7 

DX CONTAINS XF 



255 

7 






256 

0^7 

; COMPUTE X 

INTERPOLATION FACTOR 


0113 A 10600 

R 

i-J/ 

258 

7 

MOV 

AXvXSAV 

* 

* 

AX CONTAINS X LOOK-UP VALUE 

0116 2BC2 



SOB 

AX , DX 

t 

1 

SUBTRACT XF 

0118 8BDA 


260 

MOV 

BX.DX 

• 

MOVE XF TO BX REGISTER 

Oil A 8BB0 


261 

MOV 

DX,AX 

• 

> 

MOVE (XA-XF) TO DX 

011C 8BG5 


262 

MOV 

AX, BP 

■ 

? 

MOVE XG TO AX REGISTER 

01 IE 2BC3 


263 

SUB 

AX, BX 

• 

1 

SUBTRACT XF 

0120 8BC8 


264 

MOV 

CX, AX 

i 

? 

STORE RESULT IN CX 

0122 880000 


265 

MOV 

AX,0 

i 

y 

CLEAR LEAST SIG. PORTION OF DIVIDEND 

0125 B1EA 


266 

SHR 

BX, 1 

i 

* 

RIGHT SHIFT 2 TO AVOID DIVIDE OVERFLOW 

0127 D1B8 


267 

RCR 

AX, 1 



0129 D1EA 


268 

SHR 

DX, 1 



oi 2 B runs 


269 

RCR 

AX, 1 



01 2D F7F1 


270 

DIV 

CX 

» 

DIVIDE BY INTERPOLATION FACTOR 
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LOG OBJ 


LINE SOURCE 



012F A30C00 

R 

271 

MOV 

DDXiAX ; 

DDX CONTAINS X INTERPOLATION FACTOR 



272 ; 






273 : 

HXOFF CONTAINS THE HIGH X OFFSET IN THE HIGH X VECTOR (Z HIGH VALUE) 



274 ; 

LXOFF CONTAINS THE HIGH X OFFSET IN THE LOW X VECTOR (Z LOW VALUE) 



275 ? 

DBXZ CONTAINS THE Z INTERPOLATION FACTOR 



276 ! 

BDX CONTAINS THE X INTERPOLATION FACTOR 



277 ; 






278 ; 






279 ; 

COMPUTE Y HIGH LOCATIONS Y OFFSET = YBASE+X HIGH OFFSET 



280 ! 




0132 5B 


281 

POP 

BX 

BX CONTAINS THE Y TABLE BASE VALLE 

0133 SBF3 


282 

MOV 

31, BX 

SAVE Y TABLE BASE VALUE 

0135 3B3E0300 

R 

283 

MOV 

DI, HXOFF 

MOVE OFFSET TO DI REGISTER 

013“? 03DF 


284 

ADD 

BX.DI 

ADD OFFSET Y BASE 

013B 3BFB 


235 

MOV 

PI > EX 

BX CONTAINS HIGH Y LOCATION (YE) 

013D S3EF02 


236 

SUB 

DI, 2 

DI CONTAINS LOW Y LOCATION (YB) 

0140 891 El 800 

R 

237 

MOV 

YE,BT 

SAVE VALUE IN YE 

0144 893E1400 

R 

283 

■yoQ 

MOV 

YB.DI 

SAVE VALUE IN YB 



290 

COMPUTE Y 

LOW LOCATIONS Y OFFSET = YBASE + X LOW OFFSET 



291 

i 



0148 8BDE 


nnn 

MOV 

BX,SI 

BX CONTAINS THE Y TABLE BASE 

014A 3B3E0A00 

R 

293 

MOV 

DI, LXOFF 

MOVE X LOW OFFSET TO DI 

01 4E 03DF 


294 

ADD 

BX.DI 

ADD YBASE AND X LOW OFFSET 

0150 8BFB 


295 

MOV 

DI , BX 

BX CONTAINS HIGH Y LOCATION (YD) 

0152 83EF02 


296 

SUB 

01,2 

DI CONTAINS LOW Y LOCATION (YC) 

0155 893E1600 

R 

297 

OQO 

MOV 

YC.DI 

SAVE VALLE IN YC 



299 

COMPUTE YG = YD+DDXZ ( YE-YD ) 




300 

f 



0159 8B3E1800 

R 

301 

MOV 

DI.YE 

DI CONTAINS HIGH Y LOCATION YE 

015D 8B05 


302 

MOV 

AX.IDII 

AX CONTAINS HIGH Y LOCATION YD 

015F B1F8 


303 

SAR 

AX,1 

RIGHT SHIFT 1 IN CASE POINT 



304 



DELTA IS GREATER THAN 32767 

0161 8B17 


305 

MOV 

DX, [BX] 

MOVE YD INTO DX REG 

0163 D1FA 


306 

SAR 

DX,1 

RIGHT SHIFT TO AVOID SUBTRACT OVERFLOW 

0165 2BC2 


307 

SUB 

AX.DX 

SUBTRACT (YE-YD) 

0167 F72E0200 

R 

303 

IMUL 

DDXZ 

MULTIPLY BY Z INTERPOLATION FACTOR 

016B D1E0 


309 

SHL 

AX, 1 

SHIFT TO CORRECT FOR PREVIOUS 

016D B1D2 


310 

RCL 

DX, 1 

SUBTRACT OVERFLOW SHIFTS 

016F D1E0 


311 

SHL 

AX, 1 

LEFT SHIFT 1 TO CORRECT FOR DIVIDE SHIFT 

0171 D1B2 


312 

RCL 

DX,1 


0173 01E0 


313 

SHL 

AX, 1 

FULL LEFT SHIFT TO CORRECT 

0175 D1D2 


314 

RCL 

DX.l 

PREVIOUS DIVIDE SHIFT 

0177 0317 


315 

ADD 

DX, [BX] 

ADD YD VALUE TO RESULT 

0179 SBCA 


316 

MOV 

CX.DX 

CX CONTAINS YG 



317 

J 





313 

! COMPUTE YF = YC + DDXZ1YB - 

YC) 



319 




01 7B 8B3E1400 

R 

320 

MOV 

DI.YB 

5 DI CONTAINS HIGH Y LOCATION 

017F 8B1E1600 

R 

321 

MOV 

BX.YC 

5 BX CONTAINS LOW Y LOCATION 

0183 8B05 


322 

MOV 

AX.CDI] 

; MOVE YB TO AX REGISTER 

0185 D1F8 


323 

SAR 

AX, 1 

! RIGHT SHIFT 1 IN CASE 



324 



! POINT DELTA IS > 32767 

0187 3B17 


325 

MOV 

DX.CBXI 

; MOVE YG INTO DX REG. 


MCS-86 MACRO ASSEMBLER 

FUN3 



3 SEPT 30 PAGE 

LOC OBJ 

LINE 

SOURCE 



0189 01 FA 

326 

SAR 

DX.l 

SHIFT TO AV' 1 ID SUBTRACT OVERFLOW 

01 SB 2BC2 

327 

SUB 

AX 7 Da 

SUBTRACT YC 

018D F72E0200 R 

328 

IMUL 

odxz 

MULTIPLY BY INTERPOLATION FACTOR 

0191 D1E0 

32? 

SHL 

AX. 1 

SHIFT TO CORRECT FOR 

0193 D1D2 

330 

RCL 

DX.l 

PREVIOUS SUBTRACT SHIFTS 

0195 D1E0 

331 

SHL 

AX, 1 

SHIFT LEFT ONE TO CORRECT FOR DIVIDE SHIFT 

0197 D1D2 

OC'L 

RCL 

DX.l 

0199 01 EO 

333 

SHL 

AX, 1 

FULL LEFT SHIFT TO CORRECT 

019B D1D2 

334 

RCL 

DX.l 

PREVIOUS DIVIDE SHIFT 

0190 0317 

335 

ADD 

DX.IBXJ 

ADD YC TO RESULT 

019F 8PEA 

336 

MOV 

BP, DX 

BP CONTAINS YF 


337 

• 




338 

! COMPUTE FINAL YV = YF + DDXCYG - YF) 


339 

• 



03A1 SBC1 

340 

MOV 

AX.CX : 

AX CONTAINS YG 

01 AS DIFS 

341 

SAR 

AX, 1 

SHIFT TO AVOID OVERFLOW 

01 A5 01 FA 

342 

SAR 

DX, 1 : 

IF POINT DELTA >32767 

01A7 2BC2 

343 

SUB 

AX , DX 

SUBTRACT YF 

01 A9 F72E0C00 R 

344 

IMUL 

DDX ; 

MULTIPLY BY INTERPOLATION FACTOR 

01 AD D1E0 

345 

SHL 

AX, 1 

SHIFT TO CORRECT FOR PREVIOUS 

01 AF D1D2 

346 

RCL 

DX.l ; 

SUBTRACTION SHIFT 

01B1 D1E0 

347 

SHL 

AX, 1 ; 

LEFT SHIFT TO CORRECT FOR PREVIOUS 

01B3 D1D2 

343 

RCL 

DX, 1 ; 

DIVIDE OVERFLOW SHIFTS 

01B5 D1E0 

349 

SHL 

AX, 1 


01B7 D1D2 

350 

RCL 

DX, 1 


01 B? 0305 

351 

ADD 

OX, BP ; 

DX CONTAINS FINAL Y 


vSji 

; 




353 

5 RETURN TO 

MAIN PROGRAM 



354 

1 




355 

? PUSH PROGRAM COUNTER BACK ON STACK 

01 BE 3B2EOOOO R 

356 

MOV 

BP.SPC3 


01 BF 55 

357 

PUSH 

BP 


01C0 8BC2 

ocri 

OJO 

MOV 

AX. DX ; 

MOVE RESULT TO AX REGISTER FOR 


ncn 

OJ7 

Olfi 


j 

RETURN TO PL/M PROGRAM 

0102 C3 

v'OU 

361 

> 

RET 



— 

362 

CODE ENDS 




363 

END 
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XREF SYMBOL TABLE LISTING 


NAME TYPE VALUE ATTRIBUTES, XREFS 

'??SEG . SEGMENT SIZE=OOOOH PARA PUBLIC 

BEG . . L NEAR OOOEH CODE 84 861 

CFUN3 . L NEAR 0067H CODE 152 155# 

CGROUP. GROUP CODE 54# 56 

CLX3. . L NEAR 0034H CODE 170 176# 

CLXA3 . L NEAR 0031H CODE 163 174# 

CLXL3 . L NEAR OOBEH CODE 206 212# 

CLXLA3. L NEAR OOBBH CODE 199 210# 

CLZ3. . L NEAR 002AH CODE 106 112# 

CLZA3 . L NEAR 0027H CODE 98 110# 

CODE. . SEGMENT SIZE=01C3H PARA PUBLIC 'CODE' 54# 76 362 

CONST . SEGMENT SIZE=OOOOH --UNDEFINED- 55# 

DATA. . SEGMENT SIZE=001AH PARA PUBLIC 'DATA' 55# 59 73 

DDX . . V WORD OOOCH DATA 66# 271 344 

DDXZ. . V WORD 0002H DATA 61# 127 227 247 308 328 

BGROUP. GROUP DATA CONST STACK MEMORY 55# 56 

PL0P3 . L NEAR 0014H CODE 96# 99 

DLPXH3. L NEAR 006EH CODE 161# 164 

DLPXL3. L NEAR 00A8H CODE 197# 200 

FUN3. . L NEAR OOOAH CODE PUBLIC 53 85# 

HXOFF . V WORD 0008H DATA 64# 177 283 

IL0P3 . L NEAR 001DH CODE 92 103# 105 

ILPXH3. L NEAR 0077H CODE 153 167# 169 

ILPXL3. L NEAR 00B1H CODE 195 203# 205 

LXOFF . V WORD OOOAH DATA 65# 213 293 

MEMORY. SEGMENT SIZE=OOOOH — UNDEFINED-- 

NBYT. . V WORD 0004H DATA 62# 133 188 

PFUN3 . L NEAR OOOOH CODE PUBLIC 53 80# 

SPC3. . V WORD OOOOH DATA 60# 80 85 356 

STACK . SEGMENT SIZE=OOOOH —UNDEFINED— 

XB. . . V WORD OOOEH DATA 67# 182 239 

XC. . . V WORD 0010H DATA 68# 216 240 

XE. . . V WORD 0012H DATA 69# 181 220 

XSAV. . V WORD 0006H DATA 63# 155 192 258 

YB. . . V WORD 0014H DATA 70# 283 320 

YC. . . V WORD 0016H DATA 71# 297 321 

YE. . . V WORD 0018H DATA 72# 287 301 


ASSEMBLY COMPLETE, NO ERRORS FOUND 
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ISIS-II MCS-86 MACRO ASSEMBLER 02.1 ASSEMBLY OF MODULE NEFG 
OBJECT MODULE PLACED IN : F 1 : PNEFG . OB J 

ASSEMBLER INVOKED BY: ASM86 : FI: PNEFG. SRC XREF DATE 13 SEPT 80) 


LOG OBJ 


LINE SOURCE 


1 


3 


7 

O 

K.' 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
1? 
20 
21 
22 

23 

24 

25 

26 

27 

00 

i-'-l 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


SUBROUTINE NEFG FOR UNIVARIATE 

FUNCTIONS LOOK-UP OF NON-EQUAL INCREMENT FUNCTIONS 

USING THE SLOPE INTERCEPT METHOD FOR ADDITIONAL SPEED 

CALL SEQUENCES NEFG 

BP REGISTER' CONTAINS THE B VALLES' STARTING LOCATION 
SI REGISTER CONTAINS THE NUMBER CE 
POINTS IN THE TABLE 

DI REGISTER CONTAINS THE X TABLE BASE POINTER 
BX REGISTER CONTAINS THE X TABLE FLOATING POINTER ADDRESS 
AX REGISTER CONTAINS THE X LOOK-UP VALUE 
CALL NEFG 


THE Y HIGH VALUE OF THE CURVE IS RETURNED IN THE DX REGISTER 
THE Y LOW VALUE OF THE CURVE IS RETURNED IN THE BX REGISTER 
THE FLOATING X TABLE BASE POINTER IS AUTOMATICALLY UPDATED 




CALLNG SEQUENCE FOR F'L/M-86 CALL 
CALL(BPI,N,XLB,XFI,X) 

X = VALLE OF FUNCTION 

XFI= X TABLE FLOATING POINTER ADDRESS 

XLB= X TABLE BASE POINTER 

N = NUMBER OF POINTS IN TABLE 

BPI= STARTING LOCATION OF B VALUES 


»»»»»»»»»»»»»»BEWARE<«««««««««<««:«««:<«« 


THE X TABLE INDEX LOCATION IS AUTOMATICALLY UPDATED BY 
THE ROUTINE AND THEREFORE DOES NOT FOLLOW THE PL/M-86 
BLOCK STRICTURE RULES 


; »»»»»»»»»»»>»»»»<c««<« c<c<«:«««««««c«««c« 

: WHEN USING THE PL/M-86 A CALL MUST IMMEDIATELY BE MADE TO 

! THE SHFT SUBROUTINE TO GET THE CORRECT ANSWER IN THE FOLLOWING 
; MANNER: 


RSLT = SHFT (NO. 5 

WHERE NO. IS THE NUMBER OF FULL LEFT ARITHMETIC SHIFTS 
THAT MUST BE MADE 
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LOC OBJ 

LINE SOURCE 

) 


51 



52 



53 



54 

NOTE: DATA STORAGE FORMAT 


55 

X VALUES— (-32763 TO 32767) 


56 

M VALUES — ( -32763 TO 32767) 


57 

B VALLES- (DOUBLE PRECISION VALLES 


58 

LEAST SIGNIFICANT 2 BYTES 


5? 

MOST SIGNIFICANT 2 BYTES) 


60 i 



63 



62 i 



63 1 



64 



65 

NAME NEFG 


A A 

PUBLIC PNEFG.NEFG 


67 CGROUP GROUP CODE 


68 

A9 

ASSUME CS: CGROUP 


70 



71 

CODE SEGMENT FOR ROUTINE 


72 


— 

73 CODE SEGMENT PUBLIC 'CODE 7 


74 

PROGRAM ENTRY IF ENTERED BY PL/M-36 CALL ROUTINE 

0000 5A 

76 PNEFGs POP DX ; TEMP SAVE OF PC 

0001 58 

77 

POP AX ; X VALUE 

0002 5B 

73 

POP BX ; X TABLE FLOATING POINTER ADDRESS 

0003 5F 

79 

POP DI 5 X TABLE BASE VALLE 

0004 5E 

30 

POP SI ; NUMBER OF POINTS IN TABLE 

0005 5D 

81 

POP BP ; STARTING LOCATION OF B VALLES 

0006 52 

32 

00 

PUSH DX ; RESTORE PC TO THE STACK 


uo 

84 

PROGRAM ENTRY IF CALLED BY ASSEMBLY 


85 

C’i 

LANGUAGE ROUTINE 

0007 B90200 

oo 

37 1 

1EFG: MOV CX,2 ? PUT VALLE OF 2 IN CX REGISTER FOR INCREMENTING 

OOOA 3B17 

38 

MOV DX.CBX] ; MOVE POINTER VALUE TO DX REGISTER 

OOOC 87DA 

89 

XCHG BX»DX 5 SAVE POINTER LOC: IN DX REGISTER 

OOOE 03F1 

90 

ADD SI.CX ? ADD 2 TO THE NUMBER OF X POINTS 


91 


0010 3BFF7F 

92 

CMP AX, 32767 5 CHECK IF X VALLE IS EQUAL TO FULL SCALE 

0013 7501 

93 

JNZ CNEFG ; JUMP TO CONTINUE FUNCTION LOOK-UP 

0015 48 

94 

DEC AX ; DECREMENT BY ONE BIT TO AVOID POINTER 


95 

BEING LOST 


96 

BEGIN COMPARISON OF TABLE VALLES 


97 

» 

0016 3B01 

98 

CNEFGrCMP AX,[BX+DII ? COMPARE X VALUE WITH TABLE LOCATION 

0018 7008 

99 

JGE ILOP ; JUMP TO INCREMENT LOOP 


100 

DECREMENT LOOP 

001 A 2BD9 

101 

DLOP: SUB BX,CX 5 SUBTRACT 2 FROM X TABLE POINTER 

001C 3B01 

102 

CMP AX>[BX+BI] ? COMPARE X VALLE WITH TABLE LOCATION 

00 IE 7D0B 

103 

JGE INTRA ; JUMP TO INTERPOLATION ROUTINE 

0020 EBF8 

104 

JMP DLOP ; JUMP TO BEGINNING OF DECREMENT LOOP 


105 

INCREMENT LOOP 
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LOG OBJ LINE SOURCE 


0022 030? 

106 

HOP: ODD 

BX.CX ; 

ADD 2 TO TABLE POINTER 

0024 3R01 

107 

CMP 

AX. CBX+DI ] ; 

COMPARE X VALLE WITH TABLE LOCATION 

0026 7DFA 

103 

JOE 

I LOP ? 

JUMP TO BEGINNING OF INCREMENT LOOP 

0028 EB0390 

109 

JMP 

INTR ; 

JUMP TO INTERPOLATION ROUTINE 


110 

• 




111 

112 

113 

114 

5 BX CONTAINS THE LOCATION OF THE X TABLE VALUE 

; WHICH IS HIGHER THAN THE X LOOK-UP VALUE 

» 


? 

? COMPUTE LOCATION OF THE M SLOPE VALUE 


115 

! M LOCATION = ((X POINTER -(2 * #X POINTS) )-2 


116 

* 



002B 0309 

117 

INTR6; ADD 

BX.CX ? 

ADD 2 TO TABLE POINTER FOR EQUIVALENCE . 


118 



IF COMMING FROM DECREMENT LOOP 

002D 87DA 

119 

INTR: XCHG 

BX . DX ; 

PUT ADDRESS OF POINTER IN BX REGISTER 

002F 3917 

120 

MOV 

IBXI.DX 5 

SAVE CURRENT VALLE OF INDEX POINTER 

0031 0307 

121 

ADD 

DX.DI ; 

ADD X INDEX AND BASE VALLE 

0033 SBDA 

122 

MOV 

BX, BX ; 

MOVE TOTAL POINTER TO BX REG 

0035 DIES 

123 

SAL 

SI.l 5 

MULTIPLY NO. OF X POINTS BY 2 

0037 03DE 

124 

ABB 

bx.si ; 

ADD ABOVE TO CURRENT X TABLE POINTER 

0039 2BD? 

125 

SUB 

BX.CX ; 

DECREMENT POINTER BY TWO 


126 

; BX CONTAINS THE ADDRESS OF THE SLOPE M 

003B F72F 

127 

IMUL 

WORD PTR CBX] 

? MULTIPLY CURRENT X VALLE BY THE SLOPE 

0030 D1E0 

1 nr, 
iiO 

SHL 

AX, 1 

; FULL LEFT ARTIHMETIC SHIFT TO 

003F D1D2 

12? 

130 

RCL 

* 

DX, 1 

i BRING MULTIPLY SCALING INTO LINE 


131 

132 

133 

? AX AND DX CONTAIN THE VALUE OF SLOPE TIMES VALUE — H»X 


7 

: COMPUTE THE LOCATION OF THE B'S 


134 

: USING THE RELATION 

(M ADBR - X BASE - X BYTES) *2 + <B BASE) 

0041 2BDF 

135 

SUB 

BX.DI 

? SUBTRACT X BASE FROM M ADDR 

0043 2BDE 

136 

SUB 

BX.SI 

; SUBTRACT X « OF BYTES 

0045 D1E3 

137 

SAL 

BX, 1 

; MULTIPLY RESULT BY 2 

0047 03DD 

100 

ABB 

BX.BP 

; ADD BP TO GIVE ADDRESS OF MSB B VALLE 

0049 SBFB 

13? 

MOV 

DI.BX 

: MOVE VALUE TO DI REGISTER 

0048 03F9 

140 

ADD 

DI.CX 

i DECREMENT BY 2 AND DI CONTAINS ADDRESS OF 


141 



: LSB B VALUE 


142 

5 NOW FOR DOUBLE PRECI 

SION ADD 

004D 0305 

143 

ADD 

AX.CDI] 

? ADD LSB F VALUE TO RESULT 

004F 1317 

144 

ADC 

DX, CBX ] 

? ADD MSB B VALUE TO RESULT 

0051 8BD8 

145 

MOV 

BX, AX 

: MOVE CX TO AX TO ALLOW FOR THE 


146 



? SHFT FUNCTION CALL IN PLM PROGRAM 

0053 C3 

147 

RET 



— 

148 

CODE ENDS 




149 

END 
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MCS-36 MACRO ASSEMBLER NEFG 
XREF SYMBOL TABLE LISTING 


NAME 

TYPE 

VALUE 

ATTRIBUTES, XREFS 

■'■•SEG . 

SEGMENT 


SIZE=0000H PARA PUBLIC 

CGROUP. 

GROUP 


CODE 67# 63 

CNEFG . 

L NEAR 

001 AH 

CODE 93 93# 

CODE, . 

SEGMENT 


SIZE=0054H PARA PUBLIC 'CODE 

DLOP. . 

L NEAR 

001 AH 

CODE 101# 104 

I LOP. . 

L NEAR 

0022H 

CODE 99 106# 103 

INTR. . 

L NEAR 

002DH 

CODE 109 119# 

INTRA . 

L NEAR 

002BH 

CODE 103 117# 

NEFG. . 

L NEAR 

0007H 

CODE PUBLIC 66 37# 

F’NEFG . 

L NEAR 

OOOOH 

CODE PUBLIC 66 76# 


ASSEMBLY COMPLETE, NO ERRORS FOUND 
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JGIS-II MCS-36 MACRO ASSEMBLER V2.1 ASSEMBLY OF MODULE SHFT 
OBJECT MODULE PLACED IN : FI U, 'SHFT, OBJ 

ASSEMBLER INVOKED BY: ASMS6 :F1:NSHFT.SRC XREF DATE (3 SEFT 30) 


IOC OBJ LINE SOURCE 



1 












-3 

4 




5 

7 

THIS PROCEDURE IS USED TO 'SHIFT THE BX AND DX DOUBLE 


L 

-t 

! 

PRECISION ACCUMULATORS A SPECIFIED NNUMBER OF TIMES 


/ 

7 

THE PROCEDURE IS USED AFTER CALLING THE NEFG FUNCTION 


9 

7 

ROUTINE FROM A F'L/M-36 PROGRAM 


10 

7 



11 

7 



12 

1*1 

7 

CALL SEQUENCE: 


14 

7 

• 

7 

RSLT = SHFT (NO. ) 


15 

7 

WHERE NO. IS THE NUMBER OF FULL LEFT 


16 

t 

7 

ARITHMETIC: SHIFTS THAT MUST BE MADE 


37 

7 



33 

; 



1? 



20 

7 



21 


NAME SHFT 




PUBLIC SHFT 


23 

CGROUP 

GROUP CODE 


24 

95 


ASSUME CS: CGROUP 


26 

7 

7 


— 

27 

CODE 

SEGMENT PUBLIC -'CODE' 


28 

7 


0000 3BEC 

2? 

SHFT: 

MOV BP.SP 

0002 SB4E02 

30 


MOV CX, [PPI+2 

0005 D1E3 

31 

NXT: 

SHL BX>1 

0007 D1D2 

32 


RCL DX, 1 

0009 E2FA 

33 


LOOP NXT 

OOOB SBC2 

34 


MOV AX,DX 

GOOD C20200 

35 


RET 02H 

— 

36 

CODE 

ENDS 


37 . 


END 


PAGE 1 
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XREF SYMBOL TABLE LISTING 


NAME 


TYPE 


VALUE ATTRIBUTES, XREFS 


? ? SEG . SEGMENT 
CGRDUP, GROUP 
CODE. . SEGMENT 
NXT . . L NEAR 0005H 

SHFT. . L NEAR OOOOH 


SIZE=OOOOH PAF:A PUBLIC 
CODE 233 24 

SIZE=OOIOH PARA PUBLIC 'CODE" 

CODE 31# 33 

CODE PUBLIC 22 29# 


23# 


36 


ASSEMBLY COMPLETE. NO ERRORS FOUND 


INTELSCL, 07/29/80 08:30:46 


PAGE 002 


100 

200 C 
30 0 C 
4 0 0 C 
500 C 
6 0 0 C 
■ 70 0 C 
800 C 
900 C 

1 0 0 0 c 
1100 c 

1200 C 
130 0 C 
1400 C 
1500 C 
1600 C 
1700 C 
1800 C 
190 0 C 
2000 C 
2100 C 
2200 C 
2300 C 
2400 C 
2500 C 
2600 C 
27 00 C 
2800 C 
2900 C 
3000 C 
3100 C 
3200 C 
3300 C 
3400 C 
3500 C 
3600 C 
3700 C 
3800 C 
3900 C 
4000 C 
4100 C 
4200 C 
4300 C 
4400 
4500 
460 0 
47 00 
4800 
4900 
50 00 
5100 
5200 
5300 


SUBROUTINE INTSCL ( INFO , TITLE, XNUM , XDNOM , YNUM, YDNOM, ZNUM, ZDNOM, FUNC, FRMT, ORG) 

SUBROUTINE INTSCL ACCEPTS DATA IN ENGINEERING UNITS FOR NEFG, FUN1 , FUN1 A , FUN2 , AND FUN3 TYPE CURVES, 
COMPUTES M AND B FOR NEFG, AND SCALES ALL DATA. OUTPUT IS IN TWO FORMS: 

1) A TABULAR LIST DATASET, AND 

2) EITHER AN ABSOLUTE OR RELOCATABLE DATASET SUITABLE FOR THE 8086 ASSEMBLER. 

CALLING VARIABLES: 

INFO — ARRAY WHICH CONTAINS CIN THIS ORDER) 

1) X BORDER SPECIFICATION 

0.0 = ADD BOTH FLATS — NO EXTRAPOLATION 

1.0 = ADD LOW X FLAT ~ EXTRAPOLATE HIGH END 

2.0 = ADD HIGH X FLAT — EXTRAPOLATE LOW END 

3.0 = NO FLATS — EXTRAPOLATE BOTH ENDS 

2) X BREAKPOINT SPECIFICATION 

0.0 = SAME X BREAKPOINT (FUN1, FUN2 TYPE) 

1.0 = DIFFERENT X BREAKPOINT CFUN3 TYPE) 

3) NUMBER OF X'S PER Z CURVE 

4) NUMBER OF Z CURVES CO.O OR 1.0 = 1 Z CURVE) 

5) FIRST X VALUE FOR ALL Z CURVES 

6) SECOND X VALUE FOR ALL Z 

7) ETC. 

8) Z VALUES 

9) . FIRST Y VALUE FOR ALL Z CURVES 

1 0 ) ETC . 

TITLE -- NAME OF CURVE 

XNUM — X NUMERATOR VALUE FOR SCALING 

XDNOM — X DENOMINATOR VALUE FOR SCALING 

YNUM — Y NUMERATOR 

YDNOM ~ Y DENOMINATOR 

ZNUM — Z NUMERATOR 

ZDNOM -- Z DENOMINATOR 

FUNC -- INDICATES FUNCTION ROUTINE DATA BELONGS TO 

0 = NEFG 

1 = FUN1 
11= FUN1A 

2 = FUN2 

3 = FUN3 

FRMT — FORMAT CHOICE FOR OUTPUT DATASET 

0 = RELOCATABLE 

1 = ABSOLUTE 

ORG ~ HEXADECIMAL STARTING ADDRESS OF DATA FOR ABSOLUTE FORMAT 
DIMENSION XVALC100, 1 0 0 ) , ZVAL ( 1 0 0 ) ,YVAL(100, 1 0 0 ) , MVAL ( 1 0 0 ) ,BVAL(10Q) 

DIMENSION XSCLDC100,100),ZSCLD(100),YSCLD(100,100),MSCLD(100) , BMSB ( 1 0 0 ), BLSB ( 1 0 0 ) 

DIMENSION TITLEC2) , INFO C 500 ) , ZNAMEC 2 ) , XNAMEC 2) , YNAMEC 2 ) , Y1NAMEC 2 ) , MNAMEC 2 ) , BNAMEC 2 ) 

INTEGER XNO ,ZNO , XSUB , ZSUB , YSUB , ZSCL D, XSCLD , YSCLD , BMSB , BLSB 

INTEGER FUNC, ORG, FRMT, XLNS , B LNS , YLNS , ZLNS , STMTNO 

INTEGER TITLE,SHFT,Y1,Y2, EXTRAP, SFTR, WARN 

REAL MBIG,MSMAL , MVAL ,MSF 

REAL LSB , INFO ,MNAME 

WARN=0 


oj 

C D 


C 
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5400 C FILL NAME ARRAYS FOR OUTPUT DATASETS 
5500 C 

5600 DATA DAT AZ, DAT AX. DAT AY , DAT AY1.DAT AM »DATAB, BLANK/’ Z 

5700 C 

5800 IF (FUNC.NE.2.AND.FUNC.NE.3) GO TO 80 

5900 CALL F4MVC(DATAZ, 1 ,ZNAME, 1 » 1 ) 

6000 CALL F4MVC( TITLE, 1, ZNAME. 2. 6) 

6100 CALL F4MVC(BLANK,1, ZNAME, 8,1) 

6200 C 

6300 80 CONTINUE 

6400 CALL F4MVCC DATAX, 1 , XNAME, 1 , 1 ) 

6500 CALL F4MVC< TITLE, 1, XNAME, 2, 6) 

6600 CALL F4MVC(BLANK,1, XNAME, 8,1) 

6700 C 

6800 CALL F4MVC(DATAY, 1 , YNAME, 1 , 1 ) 

6900 CALL F4MVC( TITLE, 1,YNAME,2,6) 

7000 CALL F4MVC( BLANK, 1 ,YNAME,8, 1 ) 

7100 C 

7200 IF (FUNC.NE.O) GO TO 85 

7300 CALL F4MVC( DATAM, 1 ,MNAME, 1 , 1 ) 

7400 CALL F4MVC( TITLE, 1 ,MNAME, 2 , 6 ) 

7500 CALL F4MVCCBLANK, 1 ,MNAME,8 » 1 ) 

7600 C 

7700 CALL F4MVCIDATAB, 1 , BNAME, 1 , 1 ) 

7800 CALL F4MVC( TITLE, 1, BNAME, 2, 6) 

7900 CALL F4MVC(BLANK,1, BNAME, 8,1) 

8000 C 

8100 85 CONTINUE 

8200 IF (FUNC.NE.il) GO TO 90 

8300 CALL F4MVC( DAT AY1»1,Y1NAME,1,2) 

8400 CALL F4MVC(TITLE,1,Y1NAME,3,5) 

8500 CALL F4MVC(BLANK,1,Y1NAME,8,1) 

8600 C 

8700 90 CONTINUE 

88 0 0 C 
8900 C 

9000 C BREAK UP INFO ARRAY 
9100 C 

9200 CVEND=INFO(l) 

9300 BRKPT=INF0(2) 

9400 XN0=INF0(3) 

9500 C 

9600 C DETERMINE FLAT AND EXTRAPOLATION OPTIONS 
9700 C 

9800 FLAT=CVEND 

9900 C 

10000 IF (CVEND.EQ.0.0) EXTRAP=0 

10100 IF (CVEND.EQ.l.O) EXTRAP=2 

10200 IF (CVEND.EQ.2.0) EXTRAP=3 

10300 IF (CVEND.EQ.3.0) EXTRAP=1 

10400 ZN0=INF0(4) 

10500 C 

10600 C TEST FOR Z VALUES 


I 
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10700 
10800 
10 900 
11000 
11100 
11200 . 
11300 

moo 

1150 0 
11600 
11700 
11800 
11900 
12000 
12100 
12200 
12300 
12400 
12500 
12600 
12700 
12800 
12900 
13000 
13100 
13200 
13300 
1340 0 
13500 
13600 
13700 
13800 
13900 
14000 
14100 
14200 
14300 
14400 
14500 
14600 
14700 
14800 
14900 
150 00 
15100 
15200 
15300 
15400 
15500 
15600 
15700 
15800 
15900 


C 

IF (ZNO.EQ.O) ZNO = l 
JJ=ZNO 

IF ( BRKPT . EQ . 0 . 0 ) JJ=1 
C 

C GET X VALUES INTO ORDER BY CURVE 
C 

11 = 0 

DO 100 1=1, XNQ 
DO 100 J=1,JJ 
11 = 11+1 
XSUB=4+I I 

XVAL(J,I)=INFQ(XSUB) 

100 CONTINUE 
C 

C FILL ZVAL ARRAY IF NECESSARY 
C 

IF (ZNO.EQ.l) GO TO 300 
DO 200 1=1, ZNO 
ZSUB=5+XN0*JJ+(I-1) 

ZVAL(I)=INFO(ZSUB) 

200 CONTINUE 
C 

IF (FUNC.NE.2) GO TO 300 
C 

C SHIFT Z VALUES AND INSERT LOW AND HIGH ENDPOINTS 
C 

I I =ZNO+l 
III =ZNO 

DO 210 1=1, ZNO 
ZVAL (II ) =ZVAL( III ) 

ir=n.-i 

iii=iii-i 

210 CONTINUE 
C 

ZVAL ( 1 ) =“32768*ZNUM/ZDN0M 
JZN0=ZN0+2 

ZVAL ( JZNO ) =32767*ZNUM/ZDN0M 
C 

c 

C FILL Y ARRAY , 

C 

300 IF (ZNO.EQ.l) GO TO 310 
ZZNG=ZNO 
GO TO 320 
310 ZZNO=0 
C 

320 11=0 

DO 350 1=1, XNO 
DO 350 J=1 , ZNO 
YSUB=5+ (XNO* J J+ZZNO ) +11 
YVAL(J,I)=INFO(YSUB) 

11 = 11+1 
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16000 
16100 
16200 
16300 
16400 
16500 
16600 
16700 
16800 
16 9 0 0 
170 00 
17100 
17200 
17300 
17400 
17500 
176 00 
17700 
17800 
17900 
18000 
18100 
18200 
18300 
18400 
18500 
18600 
18700 
18800 
18900 
19000 
19100 
19200 
19300 
19400 
19500 
196 00 
19700 
19800 
1 9900 
20000 
20100 
20200 
20300 
20400 
20500 
“ 20600 
20700 
20800 
20900 
21000 
21100 
21200 


350 CONTINUE 
C 
C 

C EXTRAPOLATE AND ADD FLATS IF NECESSARY 

C 

C 

C TEST FOR TYPE OF EXTRAPOLATION AND PERFORM IT 
C 

IF (EXTRAP. EQ.Q) GO TO 3000 
IF ( EXTRAP . EQ . 2 ) GO TO 2000 
C 

C EXTRAPOLATE LOW END OF CURVES 

C 

C 

C SHIFT X VALUES AND INSERT LOW VALUE 
C 

1000 DO 1110 J=l,ZNO 
II=XNO+l 
III =XNO 

DO 1100 1=1, XNO 
XVALC J , II )=XVAL( J, III) 

11=11-1 
111=111-1 
1100 CONTINUE 

XVAL(J,1) =-32768#XNUM/XDNOM 
1110 CONTINUE 
C 

C SHIFT Y VECTOR VALUES 
C 

DO 1210 J=1 ,ZNO 
I I =XNO+l 
I I I =XNO 

DO 1200 1=1, XNO 
YVAL(J,II) =YVAL ( J , III ) 

11=11-1 
111=111-1 
1200 CONTINUE 
C 

C COMPUTE EXTRAPOLATED Y VALUE AND PUT IN ARRAY 

° YLOEXT = - ( ( (XVALC J, 3) -XVALC J, 1) ) * ( YVAL ( J , 3 )-YVAL C J > 2) )/(XVAL ( J, 3) -XVALC J ,2) ) )-YVAL(J,3) ) 

YVAL C J , 1 ) =YLOEXT 
1210 CONTINUE 
C 

C INCREMENT XNO TO REFLECT CHANGE 
C 

XNO=XNO+1 

C 

IF C EXTRAP . EQ . 3 ) GO TO 3000 
C 

C EXTRAPOLATE HIGH END OF CURVES 

C 

C 


CO 
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21300 C ADD X VALUE TO HIGH END 
21400 C 

21500 2000 XNO=XNO+l 

21600 C 

21700 DO 2010 J=l,ZNO 

21800 XVAL(J,XNO)=32767*XNUM/XDNOM 

21900 2010 CONTINUE 

22000 C 

22100 C COMPUTE Y HIGH EXTRAPOLATION AND PUT IN ARRAY 
22200 C 

22300 I=XNO-l 

22400 I I -XN0-2 

22500 DO 2020 J=l,ZNO 

22600 YHIEXT=(XVAL(J,XNO)-XVAL(J,II))*(YVAL(J,I)-YVAL(J,II))/(XVAL(J,I)-XVAL(J,II) )+YVALCJ,II) 

22700 YVAL(J,XNQ)=YHIEXT 

22800 2020 CONTINUE 

22900 C 

23000 C IF BOTH ENDS EXTRAPOLATED, SKIP FLATS. 

23100 C 

23200 IF (EXTRAP. EQ.l) GO TO 4500 

23300 C 

23400 3000 IF ( FL AT . EQ . 3 . 0 ) GO TO 4500 

23500 IF (FLAT. EQ. 2.0) GO TO 4000 

23600 C 

23700 C SHIFT X AND Y VALUES AND ADD FLATS 
23800 C 

23900 DO 3020 J=1,ZNQ ^ 

24000 II=XNO+l w 

24100 III =XNO 

24200 DO 3010 1=1, XNO - 

24300 XVAL(J, II) =XVAL ( J , I I I ) 

24400 YVAL(J,II)=YVAL(J, III) 

24500 11=11-1 

24600 111=111-1 

24700 3010 CONTINUE 

24800 XVAL(J,1)=-32768*XNUM/XDN0M 

24900 YVAL ( J , 1 ) =YVAL ( J , 2 ) 

25000 3020 CONTINUE 

• 25100 C 

25200 C INCREMENT XNO 
25300 C 

25400 XNO=XNO+l 

25500 C 

25600 IF (FLAT. EQ. 1.0) GO TO 4500 

25700 C 

_ 25800 C ADD UPPER FLATS TO X AND Y CURVES 

25900 C 

26000 4000 XNO=XNO+l 

26100 C 

26200 I=XNO-l 

26300 DO 4010 J=l,ZNO 

26400 XVAL(J,XN0)=32767*XNUM/XDN0M 

26500 YVAL ( J ,XNO ) =YVAL ( J , I ) 
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26600 
26700 
26800 
26900 
27000 
27100 
27200 
27300 
27400 
27500 
27 6 0 0 
27700 
27800 
27900 
28000 
28100 
28200 
28300 
28400 
28500 
28600 
28700 
28800 
28900 
29000 
29100 
29200 
; 29300 
29400 
29500 
296 00 
29700 
29800 
29900 
30 0 0 0 
30100 
30200 
30300 
30400 
30500 
30600 
30700 
30800 
30900 
- 31000 

31100 
31200 
31300 
31400 
31500 
31600 
31700 
31800 


4010 CONTINUE 
C 

4500 IF (FUNC.NE.2) GO TO 5000 
C 

C SHIFT Y VALUES AND ADD ENDPOINTS 
C 

ZN0=ZN0+2 

C 

JZN0=ZN0-2 

JJ=ZN0-1 

J J J= JZNO 

DO 370 J=I 9 JZNO 

DO 360 1=1, XNO 

YVAL(JJ,I)=YVAL(JJJ,I) 

360 CONTINUE 
JJ=JJ-1 
JJJ=JJJ-1 
370 CONTINUE 
C 

DO 380 1=1, XNO 
YVALClfI)=YVAL(2,I) 

380 CONTINUE 
C 

J=ZNO-l 

DO 390 1=1, XNO 
YVAL(ZNG,I)=YVAL( J,I) 

390 CONTINUE 
C 
C 

C DONE WITH EXTRAPOLATION AND FLATS. 

C 

C******************************************************* 

C 

C SCALING SECTION 
C 

C******************************************************* 

c 

c COMPUTE SCALE FACTORS AND MAX/MIN VALUES 
C 

5000 CONTINUE 
C 

XSF=XDNOM/XNUM 
YSF=YDNOM/YNUM 
IF (ZNUM.EQ.0.0) ZNUM=1.0 
ZSF=ZDNOM/ZNUM 
C 

XMAX=32767*XSF 

XMIN=-32768#XSF 

YMAX=32767*YSF 

YMIN=-32768*YSF 

ZMAX=32767*ZSF 

ZMIN=-32768*ZSF 

C 
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31900 
32000 
32100 
32200 
32300 
32400 
32500 
32600 
32700 
32800 
32900 
330 00 
33100 
33200 
33300 
33400 
33500 
33600 
337 00 
33800 
33900 
340 0 0 
34100 
34200 
34300 
34400 
34500 
346 00 
34700 
34300 
34900 
350 00 
35100 
35200 
35300 
35400 
35500 
35600 
35700 
35800 
35900 
36 000 
36100 
36200 
36300 
36400 
36500 
36600 
367 00 
36800 
36900 
370 00 
37100 


C SCALE Z VALUES IF NECESSARY 
C 

IF (ZNO.EQ.l) GO TO 5100 
C 

DO 5010 J = 1 > ZNO 
C 

C ROUND Z VALUES 

IF (ZVAL(J).LT.O.O) GO TO 5001 
ZSCLDC J ) =ZVAL ( J ) *ZSF+0 .5 
GO TO 5005 

5001 ZSCLD(J)=ZVAL(J)*ZSF-0.5 


C 

5005 CONTINUE 

IF (ZSCLDC J) . LE.ZMAX. OR. ZSCLDC J) .GE.ZMIN) GO TO 5006 
IF (WARN . EQ. 0 ) WRITE(10,1) 

1 FORMAT ( 1H1 ) 

WARN=1 

WR ITEC10, 5) ZVAL(J),ZSCLD(J), (TITLECN), N=l>2) 

5 FORMAT C IX , ****** THE Z VALUE *,F13.7,' AS SCALED TO ’>16 
* CURVE 1 ,2A4, ’ . ****** ) 

C 

5006 IF CZSCLDCJ) .LE. 32767) GO TO 5007 
IF (WARN . EQ . 0 ) WRITEC10,!) 

WARN=1 

WRITE(10,6)ZSCLD( J) , (TITLECN) ,N=1,2) 

6: FORMATC1X, ****** WARNING: THE Z VALUE *,I6,' FOR CURVE ' 
ZSCLDC J ) =32767 
C 

5007 IF (ZSCLDC J) .GE. -32768) GO TO 5010 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,7)ZSCLD( J) , (TITLECN) ,N=1,2) 

7 F0RMAT(1X> ****** WARNING: THE Z VALUE *,I6,* FOR CURVE * 
ZSCLDC J)=-32768 
C 

5010 CONTINUE 
C 

C SCALE X AND Y VALUES 
C 

5100 DO 5200 J=l,ZNO 
DO 5200 1=1, XNO 


, * IS OUT OF RANGE IN- 


,2A4,* HAS BEEN SET TO 32767.*****') 


, 2A4 , * HAS BEEN SET TO -32768.*****') 


C 

C ROUND AND SCALE X VALUES 

IF (XVAL(J,I).LT.0.0) GO TO 5101 
XSCLDC J,I)=XVAL( J, I)*XSF+0 .5 
GO TO 5102 

5101 XSCLD(J,I)=XVAL(J, I)*XSF-0.5 
C 

C ROUND AND SCALE Y VALUES 

5102 IF (YVALCJ,I).LT.O.Q) GO TO 5103 
YSCLDC J,I)=YVAL(J, I)*YSF+0 .5 

GO TO 5105 

5103 YSCLDC J, I )=YVALCJ, I)*YSF-0.5 


cn 
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37200 
37300 
37400 
37500 
37600 
37700 
37800 
37900 
380 00 
38100 
38200 
38300 
38400 
33500 
33600 
33700 
33800 
38900 
39000 
39100 
39200 
39300 
39400 
39500 
39600 
39700 
39800 
39900 
40000 
40100 
40200 
40300 
40400 
40500 
406 00 
40700 
40800 
40 900 
41000 
41100 
41200 
41300 
41400 
41500 
41600 
41700 
41800 
41900 
42000 
42100 
42200 
42300 
42400 


C 


5105 CONTINUE 

IF (XSCLD(J,I).LE.XMAX.OR.XSCLD(J,I).GE.XMIN) GO TO 5110 
IF ( WARN . EQ . 0 ) WRITEdQ,!) 

WARN=1 

WRITE(10,10)XVAL(J,I),XSCLD(J,I),(TITLE(N) ,N=1,2) 

10 FORMATdX, 'xxxxx THE X VALUE ',F13.7, f AS SCALED TO f ,I6, 
K CURVE f 2A4 f xx x x x 1 ) 

5110 IF (YSCLDCJiI)! LE. YMAX . OR . YSCLDC J # I) • GE. YMIN) GO TO 5120 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,15)YVAL( J,I)#YSCLD( J,I) , (TITLE(N) ,N=1 , 2) 

15 FORMATdX, 'xxxxx THE Y VALUE *,F13.7, f AS SCALED TO r ,I6, 
* CURVE 1 ,2A4, 1 . *****') 

C 

5120 IF CXSCLD(J,I) .LE. 32767) GO TO 5125 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,16)XSCLD(J, I) , ( TITLE(N) , N=1 , 2) 

16 FORMATdX, '***** WARNING: THE X VALUE SI6,' FOR CURVE 
XSCLD(J,I)=32767 

C 

5125 IF (XSCLD(J,I) .GE.-32768) GO TO 5130 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,17)XSCLD( J, I), (TITLE(N) , N=1 , 2) 

17 FORMAT (IX, 'xxxxx WARNING: THE X VALUE M6,' FOR CURVE 
XSCLD(J,I)=-32768 

C 

5130 IF (YSCLD(J> I) . LE. 32767) GO TO 5135 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,18)YSCLD(J, I), (TITLE(N) ,N=1,2) 

18 FORMATdX, 'xxxxx WARNING: THE Y VALUE M6, f FOR CURVE 
YSCLD(J,I)=32767 

C 

5135 IF (YSCLD(J,I).GE.-32768) GO TO 5200 
IF (WARN . EQ . 0 ) WRITE(10,1) 

WARN=1 

WRITE(10,19)YSCLD(J, I) , (TITLE( N ) , N=1 , 2) 

19 FORMATdX, '*#**# WARNING: THE Y VALUE M6,' FOR CURVE ', 
YSCLD(J,I)=~32768 

C 

5200 CONTINUE 
C 

C COMPUTE M AND B VALUES FOR NEFG AND SCALE 
C 

IF (FUNC.NE.O) GO TO 6000 


1 IS OUT OF RANGE IN- 


* IS OUT OF RANGE IN- 


2A4 , ' HAS BEEN SET TO 32767 . ***** • ) 


2A4, T HAS BEEN SET TO -32768 . ***** f ) 


2A4,' HAS BEEN SET TO 32767 . ***** ' ) 


2A4 , 9 HAS BEEN SET TO -32768 . ***** 1 ) 


MBNO=XNO-l 
DO 5510 I=l,MBNO 
11 = 1+1 

MVAL(I)=(YVAL(1,II)-YVAL(1,I))/(XVAL(1,II)-XVAL(1,I) ) 


CD 
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4 250 0 
42600 
42700 
42800 
42900 
43000 
43100 
43200 
43300 
43400 
43500 
43600 
43700 
43800 
43900 
44000 
44100 
44200 
44300 
44400 
44500 
44600 
44700 
44800 
44900 
45000 
45100 
45200 
45300 
45400 
45500 
45600 
45700 
45800 
45900 
46000 
46100 
46200 
46300 
46400 
46500 
46600 
46700 
46800 
46900 
47000 
47100 
47200 
47300 
47400 
47500 
47600 
47700 


BVAL(I)=YVAL(1,I)-MVAL(I)*XVAL(1,I) 

5510 CONTINUE 
C 

C COMPUTE PROPER SHIFT VALUE FOR M AND B SCALE FACTORS 
C 

MBIG=-32768 

BBIG=-32768 

MSMAL=32767 

BSMAL=32767 

C 

DO 5520 1=1 , MBNO 

IF (MVAL(I).GT.MBIG) MBIG=MVAL ( I ) 

IF (BVALCI).GT.BBIG) BBIG=BVAL ( I ) 

IF (MVAL ( I ) . LT .MSMAL ) MSMAL =MVAL ( I ) 

IF (BVAL(I).LT.BSMAL) BSMAL =B VAL ( I ) 

5520 CONTINUE 
C 

MSF=32768*YSF/XSF 

BSF=YSF 

C 

DO 5535 1=1,9 
SHFT=2**(I-1) 

SFTR=I-1 

MTST1=MBIG*MSF/SHFT 

MTST2=MSMAL*MSF/SHFT 

BTST1=BBIG*BSF/SHFT 

BTST2=BSMAL*BSF/SHFT 

C 

IF (MTST1. LE. 32767. AND. MTST2.GE. -32768) GO TO 5530 
GO TO 5535 
C 

5530 IF (BTST1.LE. 32767. AND. BTST2.GE. -32768) GO TO 5536 
C 

5535 CONTINUE 
C 

C IF SHIFTING 8 DOESN'T PUT THE SCALED VALUES WITHIN RANGE, 

C OUTPUT ERROR MESSAGE AND CONTINUE 
C 

IF (WARN . EQ . 0 ) WRITEC10,1) 

WRITE(10,20)(TITLE(N),N=1,2) 

20 FORMATC1X, 'xhkkk WARNING: THE M AND B SCALE FACTORS HAVE BEEN SHIFTED 8 (DIVIDED BY 256),'/- 
XIX,’ AND ONE OR BOTH ARE STILL OUT OF RANGE IN CURVE f ,2A4,'. *****') 

C 

C 

5536 MSF=MSF/SHFT 
BSF=BSF/SHFT 

C 

C SEPARATE B INTO MSB AND LSB 
C 

DO 5545 1=1, MBNO 
C 
C 

c 


-J 
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47800 

47900 

48000 

48100 

48200 

48300 

48400 

48500 

48600 

48700 

48800 

48900 

49000 

49100 

49200 

49300 

49400 

49500 

49600 

49700 

49800 

49900 

50000 

50100 

50200 

50300 

50400 

50500 

50600 

50700 

50800 

50900 

51000 

51100 

51200 

51300 

51400 

51500 

51600 

51700 

51800 

51900 

52000 

52100 

52200 

52300 

52400 

52500 

52600 

52700 

52800 

52900 

53000 


C ROUND AND SCALE B VALUES 
C 

IF CBMSBCD.LT. 0) GO TO 5537 

BSCLD=BVALCI)*BSF 

BMSBCI)=INTCBSCLD> 

LSB=BSCLD-BMSBCI) 

GO TO 5538 

5537 BSCLD=BVALCI)XBSF 
BMSBCI)=INTCBSCLD) 

LSB-BSCLD-BMSB Cl) 

5538 BLSBCI)=LSB#65536 
C 

C 

C ROUND AND SCALE M VALUES 

IF CMVALCD.LT. 0.0) GO TO 5539 
MSCLDCD=MVALCD*MSF+0.5 
GO TO 5545 

5539 MSCLDCD=MVALCDXMSF-0.5 

c 

5545 CONTINUE 
C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

C OUTPUT SECTION 
C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C OUTPUT NEFG TABLE AND DATASET 
C 

IF C WARN . EQ . 0 ) WRITEC10,1) 

WRITEC 1 0, 5550 )C TITLECI), I =1,2) 

5550 F0RMATC/,1X, ’NEFG CURVE ’,2A4,/) 

C 

WRITEC10,5551)XSF,XDNOM,XNUM,YSF,YDNOM,YNUM,MSF,SFTR,BSF,SFTR 

5551 FORMAT C 5X, ’ X SCALE FACTOR: f ,G12.6,’ C G12 . 6 G12 . 6 ,’)’,//, - 

*5X , 1 Y SCALE FACTOR: »,G12.6,’ C ’, G12 . 6 G12 . 6 ,’)’,//, - 

*5X,’M SCALE FACTOR: ' , G1 2 . 6 , 5X , 1 M SHIFT FACTOR: Ml,//,- 

*5X,’B SCALE FACTOR: ' , G12 . 6 , 5X, ’ B SHIFT FACTOR: ’,11,///) 

WRITEC 10 , 5552 ) 

5552 F0RMATC14X, ’X' ,13X, ’SCALED X ’ , 13X, ’ Y ’ , 13X , ' SCALED Y’//) 

C 

DO 5553 1=1, XNO 

5553 WRITEC10,5554)XVALC1,I),XSCLDC1,I),YVALC1,D,YSCLDC1,I) 

5554 FORMAT C2CSX,G13.6, 8X, 16)) 

WRITEC 10 , 5555 ) 

5555 FORMATC//9X, ’M’ ,9X, ’SCALED M’ , 9X, ’ B 9X, ’ SCALED B CMSB )’, 3X, ’ SCALED B CLSB)’) 
DO 5556 1=1 ,MBNO 

5556 WRITEC10,5557)MVALCI),MSCLDCI),BVALCD ,BMSB CI),BLSBCI) 

5557 FORMAT C3X,G13.6,4X,I6,4X,G13.6,7X,I6,11X,I6) 

C 

C TABLE DONE. NOW FOR DATASET. 

C 

WRITEC 20, 5573 )CTITL ECI ) , 1=1,2) 


CD 
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53100 
53200 
53300 
534 0 Q 
53500 
53600 
53700 
53800 
53900 
54000 
54100 
54200 
54300 
54400 
54500 
54600 
54700 
54300 
54900 
55000 
55100 
55200 
55300 
55400 
55500 
55600 
55700 
55800 
55900 
56000 
56100 
56200 
56300 
56400 
56500 
56600 
56700 
56800 
56900 
57000 
57100 
57200 
57300 
57400 
57500 
57600 
57700 
57800 
57900 
58000 
58100 
58200 
58300 
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5573 FORMAT ( 1H1 , f ; DATA FOR NEFG CURVE ? ,2A4) 

C 

C RELOCATABLE OR ABSOLUTE? 

C 

IF CFRMT.EQ.O) GO TO 5575 
WRITE(20,5574)ORG 

5574 FORMAT ( 1H ,8X, ' ORG 1 , 5X, 14 > T H ? ) 

C 

C WRITE X 
C 

5575 CONTINUE 

CALL WRTDTA(XSCLD,XNO,XNAME,0,1) 

C 

C WRITE M 
C 

5700 MLNS=MBN0/5 

IF CMLNS*5.NE.MBN0) MLNS=MLNS+1 
IF (MBN0.LT.5) GO TO 5730 

WRITEC 20, 5705 )(MNAME(N),N=1, 2), (MSCLDC I ), 1=1 , 5) 

5705 FORMAT C 1H ,2A4, * DW' , 6X, 4C 16 , 1 , ’ ) , 16 ) 

IF (MLNS.EQ.l) GO TO 5800 
KK = 6 
KKK=10 
I I =XNO 

DO 5725 I=2,MLNS 
11=11-5 

IF (II.LT.5) GO TO 5715 & 

WRITEC 20, 5710) (MSCLDC K) ,K=KK,KKK) 

5710 FORMAT ( 1H ,8X, ’ DW ? , 6X, 4 ( 16 , 1 , 1 ) , 16 ) 

KK=KK+5 
KKK=KKK+5 
GO TO 5725 
C 

5715 III=II-1 

KKK=KK+I I I 

IF (II.EQ.4) WRITEC 20, 5720) (MSCLDC K),K=KK,KKK) 

5720 FORMAT ( 1H ,8X, ' DW T , 6X, 3 ( 16 , *, ’ ) , 16 ) 

IF CII.EQ.3) WRITE(20,5721)(MSCLD(K),K=KK,KKK) 

5721 FORMAT ( 1H ,8X , f DW f , 6X , 2 C 16 , ■, 1 ) , 16 ) 

IF (II.EQ.2) WRITE(20,5722)(M5CLDCK),K=KK,KKK) 

5722 FORMAT ( 1H ,8X, T DW ? , 6X, 16 , 1 , M6 ) 

IF CII.EQ.l) WRITE(20,5723)(MSCLD(K),K=KK,KKK) 

5723 FORMAT C 1H ,8X, 1 DW S 6X, 16 ) 

C 

5725 CONTINUE 
C 

GO TO 5800 
C 

5730 CONTINUE 

IF (MBN0.EQ.4) WRITEC 20 , 5035 ) (MNAMECN ) , N=1 , 2) , (MSCLDC I ) , 1=1 , 4 ) 

5035 FORMAT ( 1H ,2A4, ’ DW f , 6X, 3( 16 , f , 16 ) 

IF (MBNO . EQ * 3 ) WRITEC 20 , 50 36 ) CMNAMEC N ) , N = 1 , 2 ) , (MSCLDC I ) , I = 1 , 3 ) 

5036 FORMAT ( 1H ,2A4 , 1 DW f , 6X, 2C 16 , f , * ) , 16 ) 
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58400 IF CMBN0.EQ.2) WRITEC20 , 5037 ) CMNAMECN) ,N = 1 , 2) , CMSCLDC I ) , 1=1 , 2) 

58500 5037 FORMAT ( 1H , 2A4 , 1 DU 1 , 6X, 16 , 1 , 1 , 16 ) 

5S60 0 IF (MBNO.EQ.l) WRI TEC 20 , 5 0 38 )MSCLD ( 1 ) 

58700 5038 FORMAT ( 1H , 2A4, r DW f , 6X, 16 ) 

58800 C 

58900 C WRITE B 
59000 C 

59100 5800 BLNS=MBN0/3 

59200 IF (BLNSX3.NE.MBNG) BLNS=BLN5+1 

59300 IF (MBNO.LT.3) GO TO 5840 

59400 WRITE(20,5805)(BNAME(N) ,N=1 ,2) ,BMSB(1) ,BLSB( 1) ,BMSB(2) ,BLSB(2) ,BMSB(3) ,BLSBC3) 

59500 5805 FORMAT C 1H , 2A4 , 1 DW f , 6X, 5 C 16 , 1 , ' ) , 16 ) 

59600 IF (BLNS.EQ.l) GO TO 9999 

59700 KK = 4 

59800 I I =MBNO 

59900 C 

60000 DO 5835 I=2,BLNS 

60100 11=11-3 

60200 IF CII.LT.3) GO TO 5815 

60300 WRITE(20,5810)BMSB(KK) , BLSB ( KK ) , BMSB C KK+1 ) , B LSBC KK+1 ) , BMSBC KK+2 ) , BLSB C KK+2 ) 

60400 5810 FORMAT ( 1H > 8X, 1 DW’ , 6X, 5 C 16 , f ) , 16 ) 

60500 KK=KK+3 

60600 GO TO 5835 

60700 C 

60800 5815 III=II-1 

60900 IF (II.NE.l) GO TO 5825 

61000 WRITEC20 , 5820 ) BMSB (KK) » BL SBC KK) cn 

61100 5820 FORMAT ( 1H ,8X, ’ DW T , 6X, 16 , ’ , 1 , 16 ) O 

61200 GO TO 9999 

61300 C 

6140 0 5825 WRI TEC 2 0 , 5830 )BMSB (KK) , BLSBCKK) , BMSB ( KK+1 ) , BLSBCKK+l ) 

61500 5830 FORMAT ( 1H , 8X , 1 DW 1 , 6X, 3 ( 1 6 , 1 , ’ ) , 16 ) 

6160 0 C 

61700 5835 CONTINUE 

61800 GO TO 9999 

61900 C 
62000 C 

62100 5840 IF (MBNO.NE.l) GO TO 5850 

62200 WRITE(2Q,5845)CBNAME(N) , N=1 , 2 ) , BMSB ( 1 ) ,BLSB(1) 

62300 5845 FORMAT ( 1H , 2A4 , 1 DW ’ , 6X , 16 , 1 , 1 , 16 ) 

62400 GO TO 9999 

62500 C 

62600 5850 WRITEC20 , 5855) CBNAMECN) ,N=1 ,2) ,BMSB(1) ,BLSB(1) , BMSB (2) ,BLSB(2) 

62700 5855 FORMAT ( 1H , 2A4 , 3( 16 , f , ’ ) , 16 ) 

62800 GO TO 9999 

62900 C 

63000 C END NEFG TABLE AND DATASET. STOP. 

63100 C 
63200 C 

63300 C FUN1 TABLE AND DATASET 
63400 C 

63500 6000 IF (FUNC.NE.l) GO TO 7000 

63600 C 


i 
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63700 
63300 
63900 
64000 
64100 
64200 
64300 
64400 
64500 
64600 
64700 
64800 
64900 
65000 
65100 
65200 
65300 
65400 
65500 
65600 
65700 
65800 
65900 
66000 
66100 
66200 
66300 
6 6 4 0 0 
66500 
66600 
66700 
66800 
66900 
67000 
67100 
67200 
67300 
67400 
67500 
67600 
67700 
67800 
67900 
68000 
68100 
68200 
68300 
68400 
68500 
68600 
68700 
63800 
68900 


IF (WARN . EQ . 0 ) WRITECIO',1) 

WRITE(10,6005)(TITLE(I) , 1=1 , 2) ,XSF, XDNOM,XNUM, YSF, YDNOM, YNUM 
6005 FORMAT (/,1X,4X, f FUN1 CURVE T ,2A4,//~ 

*5X,'X SCALE FACTOR: »,G12.6, f (', G12 . 6 , 1 /», G12 .6 ,’)’,//, - 

*5X,'Y SCALE FACTOR: »,G12.6,' ( ■ , G12 . 6 , T /', G12 . 6 ,')',/// ) 

C 

WRITE(10,6010) 

6010 FORMAT ( 14X, T X 1 , 13X, ’SCALED X 1 , 13X, 1 Y M3X , * SCAL ED Y*//) 

DO 6015 1=1, XNO 

6015 URITE(10,602Q)XVAL<1,I),XSCLD(1,I),YVAL(1,I),YSCLDC1,I) 

6020 FORMAT (2(8X,G13.6,8X,I6)) 

C 

C TABLE DONE. NOW FOR DATASET. 

C 

WRITE(20,6025)(TITLE(I) , 1=1,2) 

6025 FORMAT ( 1H1 , ' ; DATA FOR FUN1 CURVE T ,2A4) 

C 

C RELOCATABLE OR ABSOLUTE? 

C 

IF (FRMT.EQ.O) GO TO 6035 
WRITE(2Q, 6030)0 RG 

6030 FORMAT ( 1H ,8X, ' ORG 1 , 5X, 14 , 1 H 1 ) 

C 

C WRITE X AND Y 
C 

6035 CONTINUE 

CALL WRTDTA(XSCLD,XNO,XNAME, 0,1) 

C 

CALL WRTDTA(YSCLD,XNO,YNAME,0,1) 

C 

C END FUN1 TABLE AND DATASET. STOP. 

C 

GO TO 9999 
C 

C FUN1A TABLE AND DATASET. 

C 

7000 IF (FUNC.NE.il) GO TO 8000 
C 

IF (WARN . EQ . 0 ) WRITE(10,1) 

WRITE(10,7005)(TITLE(I),I=1,2),XSF,XDNQM,XNUM, YSF, YDNOM, YNUM 
7005 FORMAT (/ , IX, 4X, 1 FUN1 A CURVE S2A4,//,- 

*5X,'X SCALE FACTOR: »,G12.6,' ( * , G12 . 6 , »/ 1 , G12 . 6 , * ) »,// , “ 

*5X,'Y SCALE FACTOR: SG12.6,' ( f , G12 . 6 , '/ 1 , G12 . 6 , 1 ) ' ,/// ) 

WRITE( 10 , 6010 ) 

DO 7010 1=1, XNO 

7010 WRITE(10,6020)XVAL(1,I),XSCLD(1,I),YVALC1,I),YSCLD(1,I) 

C 

C TABLE DONE. NOW FOR DATASET. 

C 

WRITE(2Q,7015)(TITLE(N) ,N=1,2) 

7015 FORMAT ( 1H1 , T ; DATA FOR FUN1A CURVE f ,2A4) 

C 

C RELOCATABLE OR ABSOLUTE? 


cn 
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6 9000 
69100 
69200 
69300 
6 9400 
69500 
6 9600 
6 9700 
69800 

6 9900 
70 00 0 
70100 
70200 
70300 
70400 
70500 
7060 0 
7070 0 
70800 

7 0900 
71000 
71100 
71200 
71300 
71400 
71500 
71600 
71700 
71800 
71900 
72000 
72100 
72200 
72300 
72400 
72500 
72600 
72700 
72800 
72900 
7300 0 
73100 
73200 
73300 
73400 
73500 

" 73600 

73700 
73800 
73900 
74000 
74100 
74200 


C 

IF (FRMT.EQ.Q) GO TO 7020 
WRITEC2O,603Q)ORG 
C 

7020 CONTINUE 

CALL WRTDTA(YSCLD,XNO,Y1NAME,0*1) 

C 

C END FUN1A TABLE AND DATASET. STOP. 

C 

GO TO 9999 
C 

C FUN2 TABLE AND DATASET. 

C 

8000 IF CFUNC.NE.2) GO TO 9000 
C 

DO 8047 J=1 , ZNO 
IF (J.NE.l) GO TO 8010 
IF (WARN . EQ . 0 ) WRITE<10,1> 

WRITE(1Q,80Q1)(TITLE(N) ,N=1,2) 

8001 FORMAT (/, IX* 1 FUN 2 CURVE f *2A4*//) 
WRITE(1O,8005)XSF,XDNQM,XNUM,YSF,YDNGM,YNUM,ZSF,ZDNOM,ZNUM 

8005 FORMAT ( - 

*5X,'X SCALE FACTOR: ',G12.6,' ( ' , G12 . 6 * V f , G12 . 6 , T ) ' , // * - 

*5X, r Y SCALE FACTOR: »,G12.6, f ( ' , G12 . 6 * f / ? * G12 . 6 , 1 ) f , // * - 

*5X, 1 Z SCALE FACTOR: »,G12.6, f ( f , G12 . 6 , 1 / 1 , G12 . 6 * f ) 1 , /// ) 

C 

8010 IF (J.NE.l) GO TO 8016 

WRITE(10,8015)ZVAL(J),ZSCLD<J) 

8015 FORMAT ( 14X * ? Z = 1 * G13 . 6 , 13X* 1 SCALED Z = f ,I6,///) 

GO TO 8021 

8016 WRITE(10*8020)(TITLE(N) , N = 1 * 2 ) * ZVAL ( J ) * ZSCLDC J ) 

8020 FORMAT (1H1,/,1X, 1 FUN2 CURVE 1 , 2A4 , /// , 13X, f Z = f , G13 . 6 , 13X, f SCALED Z 
C 

C X AND Y TABLES 
C 

8021 WRITE(10*8Q25) 

8025 FORMAT (8X* ’X f *8X* f SCALED X f , 8X , f Y f * 8X , ' SCAL ED Y f ) 

C 

DO 8046 1=1, XNO 

WRITEC10,8030 )XVAL(1,I) , XSCLDC 1 * I ) , YVAL ( J,I) *YSCLD(J,I) 

8030 F0RMAT(2(3X,G13.6,3X,I6)) 

C 

8046 CONTINUE 

8047 CONTINUE 
C 

C DONE WITH Z TABLE. NOW FOR Z DATASET. 

C 

C 

C RELOCATABLE OR ABSOLUTE? 

C 

WRITE(20*8099)(TITLE(N) *N=1*2) 

8099 FORMAT ( 1H1 * 1 * DATA FOR FUN2 CURVE **2A4) 

8100 IF (FRMT.EQ.O) GO TO 8102 


f ,16,///) 


on 

DO 


I 
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74300 

74400 

74500 

74600 

74700 

7 4800 

74900 

75000 

75100 

75200 

7.5300 

75400 

75500 

75600 

7 5700 

75800 

75900 

76000 

76100 

76200 

76300 

76400 

76500 

76600 

767 00 

76800 

76900 

77000 

77100 

77200 

77300 

77400 

77500 

77600 

77700 

77800 

77900 

78000 

78100 

78200 

78300 

78400 

78500 

78600 

78700 

78800 

78900 

79000 

79100 

79200 

79300 

79400 

79500 


WRITEC20,8101)ORG 

8101 FORMAT ( 1H ,8X, * ORG 1 , 5X, 14 , f H 1 ) 

C 

C Z VECTOR 
C 

8102 ZLNS=ZN0/5 

IF (ZLNS#5.NE.ZN0) ZINS=ZLNS+1 
IF (ZNO . LT . 5 ) GO TO 8125 
IF (FUNC.EQ.3) GO TO 8106 

WRITE(20,8105) (ZNAME(N) ,N = 1 ,2) , (ZSCLD(I) , 1=1,5) 

8105 FORMAT C 1H ,2A4, f DW' , 6X, 4 ( 16 , ' , T ) , 16 ) 

GO TO 8108 

8106 WRITE(20,8107)(ZNAME(N) ,N = 1,2)> CZSCLDCI) ,1 = 1,5) 

8107 FORMAT ( 1H , 2A4 , r DW * , 6X , 4 C 16 , 1 , 1 ) , 16 ) 

IF (ZLNS.EQ.l) GO TO 8199 

C 


8108 KK=6 

KKK=1G 
II =ZNO 
C 

DO 8120 1=2, ZLNS 
11=11-5 

IF (II.LT.5) GO TO 8115 
WRITE(2Q,8110)(ZSCLD(K),K=KK,KKK) 
8110 FORMAT ( 1H ,8X, ' DW T , 6X, 4 ( 16 , f , f ) , 16 ) 
KK=KK+5 
KKK=KKK+5 
GO TO 8120 
C 


8115 III=II-1 

KKK=KK+ III 

IF (II.EQ.4) WRITE(20,5720)(ZSCLD(K),K=KK,KKK) 
IF (II.EQ.3) WRITE(20,5721) (ZSCLD(K) ,K=KK,KKK) 
IF (II.EQ.2) WRITE(2Q,5722) (ZSCLDCK) ,K=KK,KKK) 
IF (II.EQ.l) WRITE(20,5723)(ZSCLD(K),K=KK,KKK) 
C 

8120 CONTINUE 

GO TO 8199 
C 


8125 CONTINUE 
C 


IF (ZN0.EQ.4) WRITE(20,8130)(ZNAME(N),N=1, 2), CZSCLDCI), 1=1, 4) 

8130 FORMAT ( 1H , 2A4 , ' DW 1 , 6X , 3 C 16 , 1 , 1 ) , 16 ) 

IF CZNO.EQ.3) WRITE(20,8131)(ZNAMECN),N=1, 2), CZSCLDCI), 1=1, 3) 

8131 FORMAT ( 1H ,2A4, T DW T ,6X,2(I6, 1 , 1 ),I6) 

IF CZN0.EQ.2) WRITE(20,8132)(ZNAME(N),N=1, 2), CZSCLDCI), 1=1, 2) 

8132 FORMAT ( 1H ,2A4, f DW T ,6X,I6 , 9 , 1 , 16 ) 

IF CZNQ.EQ.l) WRITEC20,8133)ZSCLDC1) 

8133 FORMAT ( 1H ,2A4, ’ DW f , 6X, 16 ) 

C 


8199 IF (FUNC.EQ.3) GO TO 9100 
C 

C WRITE X AND Y 


cn 



INTELSCL, 07/29/80 08:30:46 


PAGE 017 


79600 
79700 
79800 
79900 
80000 
80100 
80200 
80300 
8040 G 
80500 
80600 
80700 
80800 
80900 
81000 
81100 
81200 
81300 
81400 
81500 
81600 
81700 
81800 
81900 
82000 
82100 
82200 
82300 
82400 
82500 
82600 
82700 
82800 
82900 
83000 
83100 
83200 
83300 
83400 
83500 
83600 
33700 
83800 
83900 
84000 
84100 
84200 
84300 
84400 
84500 
84600 
84700 
84800 


C 

8200 CONTINUE 

CALL WRTDTA(XSCLD,XNO,XNAME,l,l) 

C 

CALL WRTDTA(YSCLD,XN0,YNAME,1,ZN0) 

C 

C DONE WITH FUN2 . STOP. 

C 

GO TO 9999 
C 

C NOW FOR FUN3 TABLES. 

C 

c 

9000 CONTINUE 
C 

DO 9047 J=l,ZNO 
IF (J.NE.l) GO TO 9010 
IF (WARN . EQ . 0 ) WRITE(1Q,1) 

WRITE(10,9G01)(TITLE(N) ,N=1,2) 

9001 FORMAT (/ , IX, ’ FUN3 CURVE »,2A4,//) 

WRIT EC 10, 90 05)XSF,XDNQM,XNUM,YSF,YDNOM,YNUri,ZSF,ZDNQM,ZNUM 
9005 FORMATC- 

#5X,'X SCALE FACTOR: »,G12.6, f ( f , G12 . 6 , '/ 1 , G12 . 6 , 1 ) ' ,// , - 

*5X,'Y SCALE FACTOR: »,G12.6,’ ( 1 , G12 . 6 , ' / * , G12 . 6 , f ) 1 , // , - 

*5X,'Z SCALE FACTOR: f ,G12.6, f ( f , G12 . 6 , f / G12 . 6 , f ) f ,/// ) 

C 

9010 IF (J.NE.l) GO TO 9016 

WRITE(1Q,9015)ZVAL(J) ,ZSCLD( J) 

9015 FORMAT ( 14X, ' Z = f , G13 . 6 > 13X, * SCALED Z = f ,I6///) 

GO TO 9021 

9016 WRITE(10,9020)(TITLE(N) ,N=1,2) ,ZVAL( J) ,ZSCLD( J) 

9020 F0RMAT(1H1,/,1X, 1 FUN3 CURVE f ,2A4,13X, f Z = 1 , G13 . 6 , 13X, 1 SCALED Z 
C 

C X AND Y TABLES 
C 

9021 WRITE( 10*9025) 

9025 F0RMAT(8X, 'X' ,8X, 'SCALED X 1 , 8X , 1 Y 1 , 8X, 1 SCAL ED Y'//) 

C 

DO 9046 1=1, XNO 

WRIT E(10,9030 )XVAL ( J,I) ,XSCLD(J,I) ,YVAL(J,I),YSCLD(J,I) 

9030 FORMAT (3X,G13.6,3X,I6,3X,G13.6,3X, 16) 

C 

9046 CONTINUE 

9047 CONTINUE 
C 

C DONE WITH Z TABLE; NOW FOR Z DATASET. 

C 

C Z VECTOR 
C 

WRITE(20,9050)(TITLE(N) ,N=1»2) 

9050 FORMAT ( 1H1 , 1 ; DATA FOR FUN3 CURVE ',2A4) 

GO TO 8100 
C 


' ,16///) 


ui 
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84900 C X AND Y VECTORS 

85000 C 

85100 9100 CONTINUE 

85200 CALL WRTDTA(XSCLD,XNO,XNAME, l,ZNO) 

85300 C 

85400 CALL WRTDTA(YSCLD,XNO,YNAME, l,ZNO) 

85500 C 

85600 C END OF ROUTINE 

85700 C 

85800 9999 RETURN 

85900 END 


CD 

U1 
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100 SUBROUTINE WRTDTAC FUNARY , XPTS, DNAME, FNTN, ZPTS) 

200 C 

300 C THIS SUBROUTINE WRITES DATA VECTORS FOR DATA PROCESSED IN ROUTINE INTSCL. 

400 C 

500 C CALLING SEQUENCE: 

700 C FUNARY — TWO-DIMENSIONAL ARRAY FOR X AND Y VECTORS IN NEFG, FUN1 , FUN1 A , FUN2 , AND FUN3 

800 C XPTS — NUMBER OF X POINTS IN VECTOR 

900 C DNAME — VARIABLE NAME TO BE PLACED TO START OF DATA 

1000 C FNTN — INDICATES FUNCTION TABLE IS FOR 

1100 C 0 = NEFG, FUN1 , FUN1A 

1200 C 1 = FUN2 , FUN3 

1300 C ZPTS — NUMBER OF Z POINTS IN VECTOR 

1400 C 

1500 DIMENSION FUNARYC 100 , 100 ) , DNAME(2) 

1600 INTEGER FUNARY, XPTS , FNTN , ZPTS , ZZPTS 

1700 C 

1800 ZZPTS=ZPTS 

1900 IF (FNTN.EQ.O) ZZPTS=1 

2000 C 

2100 C 

2200 C 

2300 DO 2000 J=l, ZZPTS 

2400 C 

2450 100 CONTINUE 

2500 LNS=XPTS/5 

2600 IF (LNSK5.NE.XPTS) LNS=LNS+1 

2700 IF (XPTS . LT . 5 ) GO TO 1030 

2800 C 

2900 C WRITE FIRST LINE 

3000 C 

3100 IF (J.NE.l) GO TO 1006 

3200 WRITE(20,1005)(DNAME(N),N=1,2),(FUNARY(J,I),I=1,5) 

3300 1005 FORMAT (IN ,2A4, ' DW' , 6X, 4C 16 ,','), 16 ) 

3400 GO TO 1009 

3500 C 

3600 1006 WRITE(20»1010)( FUNARY (J, I), 1=1,5) 

3700 C 

3800 1009 IF (LNS.EQ.l) GO TO 2000 

3900 KK=6 

4000 KKK=10 

4100 II=XPTS 

4200 C 

4300 C WRITE SUBSEQUENT LINES 

4400 C 

4500 DO 1025 1=2, LNS 

4600 11=11-5 

4700 IF (II.LT.5) GO TO 1015 

4800 WRITE(20, 1010 )( FUNARYC J,K),K=KK,KKK) 

4900 1010 FORMAT ( 1H ,8X, ’ DW* , 6X,4( 16 ,’,’), 16 ) 

5000 KK=KK+5 

5100 KKK=KKK+5 

5200 GO TO 1025 


9 



WRITDATA, 07/29/80 08=29=44 


PAGE 003 


5300 C 

5400 C IF LAST LINE HAS < 5 POINTS: 

5500 C 

5600 1015 III=II-1 

5700 KKK=KK+III 

5800 IF CII.EQ.4) WRITEC 20 , 1045K FUNARY( J , K) , K=KK, KKK) 

5900 IF (II.EQ.3) WRITEC 20 , 1046 ) C FUNARYC J , K) , K=KK, KKK ) 

6000 IF (II.EQ.2) WRITEC20,1047)CFUNARYCJ,K),K=KK,KKK) 

6100 IF (II.EQ.l) WRITEC20,1048)CFUNARYCJ,K),K=KK,KKK) 

6200 C 

6300 1025 CONTINUE 

6400 C 

6500 GO TO 2000 

6600 C 

6700 C IF VECTOR HAS < 5 TOTAL POINTS 
6800 C 

6900 1030 CONTINUE 

7000 IF CJ.NE.l) GO TO 1040 

7100 C 

7200 IF CXPTS.EQ.4) WRITEC 20 , 1 035 ) ( DNAMECN) , N=1 , 2 ) , ( FUNARYC J , I ) , 1=1 , 4 ) 

73C0 1035 FORMATC1H ,2A4 ,* DW* , 6X, 3C 16 ,',»), 16 ) 

7400 IF CXPTS.EQ.3) WRITEC 20 , 1 036 )( DNAMEC N ) , N = 1 , 2 ) , C FUNARYC J , I ) , 1 = 1 , 3 ) 

7500 1036 FORMAT C 1H ,2A4, »DW 6X, 2C 16 ,','), 16 ) 

7600 IF CXPTS.EQ.2) WRITEC 20 , 1 037 ) C DNAMEC N ), N=1 , 2 ), C FUNARYC J , I ), 1=1 , 2 ) 

7700 1037 FORMAT C 1H ,2A4 , • DW * , 6X, 16 , ' , 1 , 16 ) 

7800 IF CXPTS.EQ.l) WRITEC 20 , 1 038 ) C DNAMEC N) , N=1 , 2 ), FUNARYC J , 1 ) 

7900 1038 FORMAT C 1H ,2A4 , • DW * , 6X, 16 ) oi 

8000 C ^ 

8100 GO TO 2000 

8200 C 

8300 1040 CONTINUE 

8400 C 

8500 IF CXPTS.EQ.4) WRITEC20 , 1 045 ) C FUNARYC J , I ) , 1=1 , 4 ) 

8600 1045 FORMATC1H ,8X, • DW f , 6X, 3C 16 16 ) 

8700 IF CXPTS.EQ.3) WRITEC 20 , 1 046 ) C FUNARYC J , I ), 1=1 , 3) 

8800 1046 FORMAT C 1H ,8X, • DW 6X, 2C 16 ,',’), 16 ) 

8900 IF CXPTS.EQ.2) WRITEC 20 , 1 047 ) C FUNARYC J , I ) , I =1 , 2) 

9000 1047 FORMAT C 1H ,8X, 'DW' , 6X, 16 , » , 1 , 16 ) 

9100 IF CXPTS.EQ.l) WRITEC 20 , 1048 ) FUNARYC J , 1 ) 

9200 1048 FORMAT C 1H ,8X, • DW* , 6X, 16 ) 

9300 C 
9400 C 

9500 2000 CONTINUE 

9600 9999 RETURN 

9700 END 






y 



Figure 5. - Map function interpolation diagram. 


-32768 Xj x 2 x 3 x 4 32767 

yi yi y2 y3 *4 >4 

Figure 6 . - Example of FlNl data 
storage format. 


x value -3276S x 2 x 3 x 4 32768 

y value yj yi Yi Vi y4 

y value prime y^ y^ ^2 y 4 > 4 ' 

Figure 7. - Example of FUN1 and FUN1A data storage 
format. 


x value -32768 *j * 2 X 3 x 4 32768 
m value ntj m 2 m 3 m 4 m 5 

b value b 1 b 1L bfri ty> 3 i »4 b 4l b 5 b 5L 

Figure 8 . - Example of NEFG data storage format. 


32767 


z values *32768 ij 13 z 3 
x values *32768 Xj Xj *3 x 4 ^2767 

y values yjj yjj yj 2 yi3 ^14 *14 

yn *n *12 *» *w *w 

*21 *21 *22 *23 *24 *24 

y 3 l *31 *32 *33 *34 *34 

*31 *31 *32 *33 *34 *34 

Figure 9. - Example of FUN 2 data storage format. 
(Note y:: form, where i corresponds to z 
value and j corresponds to x value. 


Z 1 

z 2 

z 3 

z 4 


-32768 

X 11 

x 12 

x 13 

32767 

-32768 

X 21 

x 22 

x 23 

32767 

-32768 

X 31 

x 32 

x 33 

32767 

-32768 

X 41 

x 42 

x 43 

32767 

*11 

*11 

*12 

*13 

*13 

*21 

*21 

*22 

*23 

*23 

*31 

*31 

*32 

*33 

*33 

*41 

*41 

*42 

*43 

*44 


Figure 10. - Example of FUN 3 data storage for- 
mat. (Note that Xj; and yj: form, where 
i corresponds to z value and j corresponds 
to x value.) 


Data/function type/XB, XBRKP, #x‘s, #z's 
x values 
z values 
y values 

Figure 11. - Generalized data input. 



Data/FUN 1/0.0. a 0, 3. 0. 0. 0 

*!• x 2' x 3 
Vl> y 2- y 3 

Figure 12. - Example of FUN1 
data input for Fortran program. 


Data/FUN 2/0. 0, 0.0, .4.0, 3.0 

Xj, x 2 , Xj, X 4 

z l- z 2 > z 3 

y11.y12.y13 

y21.y22.y23 

y 31’ y 32 y 33 

y4i.y42.y43 

Figure 13. - Example of FUN2 
data input for Fortran pro- 
gram. Mote y : : form, where 
i corresponds to x value and 
j corresponds to z value. ) 


Data/FUN 3/ 0.0, 1.0, 4.0, 3.0 
xu, x 12 . x 13 

X 2i . x 22 , x a 

x 31’ *32' x 33 
x 41’ x 42' x 43 
z l- z 2' z 3 

yn- vi2- vi3 

y21> y 22’ y 23 
y 31* y 32' y 33 
y 4 i. y42- y 43 


Figure 14. - Example of FUN3 
data input for Fortran pro- 
gram. Note yi: form, where 
i corresponds lb x value and 
j corresponds to z value.) 



Subroutine 
{ INTELSCl 

Read data elements 
from calling program 


Analyze calling parameters 
for extrapolation information 
and rearrange data 


FUN2 or 
FUN3 curves 















FUN1 CURVE DP2SSH 

x scut nr.w- «i.«oo 

Y SCALE FACTO* 1 «««. 


( 32HM 
( 327(8 . 0 


/ Aio.no 

/•. 14*000 


) 

) 


-400.000 

73.4000 

110.000 

123.400 

144.000 

132.000 

180.400 

203.000 

204.000 

218.000 
223.700 
232.500 
3f*.*88 


-32748 

4020 

0748 

1*142 

11040 

12452 

14778 

14430 

14712 

17850 

18324 

10044 

32747 


*.7400008-01 
0 .740000E-01 
1.106000 
*.100350 
8.115400 
*.115540 
1. 115240 
0.115010 
0.115000 
0.111450 
4.110000 
1.102000 
0.102000 


15545 
13545 
22116 
22305 
23475 
23443 
2340 1 
23554 
23532 
22825 
22526 
20600 
20600 


(a) FUN1. 


HEFG CURVE M1F2CV 

X SCALE FACTO*: 800.000 
Y SCALE FACTO*: 2.18453 
n SCALE FACTO*'- 80.4784 
8 SCALE FACTO*: 2.18453 


C 32000.0 / 40.0000 

( 32748.0 / 15000.0 

n SHIFT factor: 0 
8 SHIFT FACTOR: 0 


•40.0400 

25.0CC0 

30.0000 

33.0000 
40.0387 


0.000000 

04.5000 

135.100 

1.000000 


-32745 

'20000 

riooo 

000 

:74? 


o.esoooo 

-2412.50 

-3370.50 

1136.00 


0.000000 

0.000000 

462.500 

1156.00 

1138.00 


1054 

2530 

2530 


SCALED B (MB) SCALED I USB) 


-3270 

-7709 

2320 


-12032 

-370*8 

43168 


(b)NEFG. 


Figure 16. - Examples of scaling program output for FIN1 and NEFG curves. 


FUH2 CURVE PT4SCH 

X SCALE FACTO*: 11.0200 

Y SCALE FACTO*: 320.000 

Z SCALE FACTO*: 800.000 

Z * -40.0400 


< 32748.0 / 400.000 
C 32000.0 / 100.000 
( 32000.0 / 40.0000 

SCALED Z « -32748 


-400 . 000 

111.000 

145.000 

100.000 
234.000 
390.088 


SCALED X 
-32748 
0003 
13517 
14302 
10140 
32747 


1 .80000 
1.60000 


2 ! 0 0000 
4.10000 




SCALED Y 
574 
574 
406 
440 
1312 
1312 


UARHING: 
».mi WARN I NG 5 
.•••a UARHING: 

• mm UARHING: 

• ml UARHING: 


THE Y VALUE 
THE Y VALUE 
THE Y VALUE 
THE Y VALUE 
THE Y VALUE 


44262 FOR CURVE H7 
53800 FOR CURVE M7 
inis FOR CURVE H7 
42386 F0« CURVE H7 
33147 FO* CURVE H7 


FUN3 CURVE H7 


X SCALE FACTO* 5 2.13335 
Y SCALE FACTO*: 727.273 
2 SCALE FACTO*: 32.0000 


( 32000.0 
( 32000.0 
( 32000.0 


HAS BEEN SET TO 32787. 

HAS BEEN SET TO 32747. 

HAS SEEM SET TO 32747. 

4FT TO 


/ 15000.0 ) 
/ 44.0000 > 
/ 1008.00 ) 


-153400 

7014.00 

8447 .00 
*05*. 00 
1535*. 5 


Z a -100.000 

SCALED X 

-32748 
14*47 
18020 
1*324 
32747 


-40.0000 

-40.0000 

- 2.00000 

4.00000 

883598 


SCALED 2 * -3200 

SCALED Y 


- 2 * 0*1 

-290*1 


-1455 

4344 

32747 


(b) FUN 3. 


Figure 17. - Example of scaling program output for FIN2 and FIN 3 curves. 



data ro* 

FIIN2 cuovi UACC8V 

ZUACCtV 

ou 

“327*8. 

-2*1. 

80*. 2317, 


DU 

*133, 

11831. 

177*1, 327*7 

XUACCftV 

DU 

-327*8. 

*2*7. 

5173. 7*80, 11013 


DU 

1*213, 

17707, 

32767 

YUACCOV 

DU 

177 3 , 

1773. 

10781. 13001, 1*7*1 


DU 

1832* , 

111*1. 

111*1 


DU 

1773, 

1773. 

10781, 13001, 1*7*1 


DU 

1132*. 

111*1. 

111*1 


cu 

8*13. 

8*13. 

10781, 13001, 1*03* 


DU 

17883. 

no** . 

110«* 


DU 

77*1. 

77*1. 

10781. 12*33. 1*771 


CU 

1****. 

188*2. 

188*2 


DU 

77*1. 

77*1. 

10781, 111*1. 13*5* 


CU 

13311. 

1819*. 

1810* 


cu 

77*1. 

77*1 , 

13781, 11*28. 127*7 


cu 

13118. 

1*137 . 

1*137 


DU 

77*1. 

7 7*1, 

10781. 10781. 10781 


:u 

11*28. 

1*8*8. 

1*8*8 


CU 

77*1 . 

77*1. 

10781. 10*01. 10401 


DU 

10*01. 

11310. 

11510 


DU 

77*1. 

77*1 , 

10781. 10*01, 10*01 


DU 

10*01. 

11310, 

11310 


Figure 18. - Example of scaling program code 
output for FUN2 curve. 
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